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L  Overview  of  Research  Program 

High  Power  Semiconductor  Switches  (HPSS)  have  received  considerable  attention 
for  the  promise  they  offer  to  pulse  power  technologies.  This  promise  is  founded  on 
the  fact  that  with  semiconducting  materials  as  the  switching  medium,  it  is  in  principle 
possible  to  obtain:  a)  high  current  densities  over  large  areas,  b)  conductivity 
modulation  over  several  orders  of  magnitude,  and  c)  high  dielectric  field  strengths. 
Optically  or  electrically  triggered  bulk  and  junction  devices  are  being  investigated  in 
a  variety  of  geometries  and  covering  a  wide  range  of  parameters.  Semiconducting 
materials  that  have  been  considered  are  Si,  GaAs,  and  very  recently,  diamond. 
Because  of  its  high  mobility  and  hold-off  voltage,  GaAs  devices  have  been  the  focus 
of  most  of  the  effort  However,  a  major  difficulty  has  surfaced  with  regards  to  the 
forward  voltage  drop  of  GaAs  devices  in  the  on-state.  That  is,  the  phenmnenon  of 
lock-on;  upon  switching,  the  device  evolves  into  a  self-sustaining  state  characterized 

by  a  constant  voltage  with  high  average  field  (z  4  kVJcmy,  and  consequently, 

intolerable  losses.  A  number  of  possible  causes  have  been  proposed  involving  either 
contact  or  bulk  properties  of  the  devices.  However,  the  physics  behind  this 
phenomenon  is  yet  to  be  understood. 

With  AFOSR  support,  our  research  activities  have  focused  on  developing 
theoretical  models  and  experimental  diagnostic  techniques  for  investigating  this  and 
other  transient  characteristics  of  these  devices.  Our  aim  has  been  to  devetop  a 
quantitative  tmderstanding  of  the  role  of  the  various  microscopic  processes,  material 
parameters,  traps,  and  space-charge  in  shaping  the  behavior  of  the  devices.  During 
this  contract  period,  we  have  developed  a  non-equQibrium  model  for  describing  the 
dynamics  of  carriers  in  multi-valley  semiconductors.  This  model  has  been  applied  to 
GaAs  devices  subjected  to  rapidly  varying  (in  time)  fields.  A  code  developed  during 
the  previous  contract  period  has  been  used  to  obtain  a  set  of  material  parameters  for 
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GaAs  by  minimizing  the  error  between  calculated  and  measured  values  for  transport 
parameters.  The  experimental  facility  to  investigate  the  dynamics  of  switching,  in 
particular  space-charge  redistribution,  has  also  been  completed  during  this  period 
In  the  next  section  a  summary  of  the  accomplishments  made  during  the  period  of  this 
contract  is  given.  Sections  III-V  are  devoted  to  a  discussion  of  the  results. 

IL  Review  of  Achievements  and  Activities  During  the  Contract  Period 

During  this  funding  period,  we  have  accomplished  the  following: 

(1)  The  concept  of  a  macrokinetic  distribution  (MKD)  developed  in  the  previous 
contract  period  to  describe  the  non-equilibrium  dynamics  of  carriers  in  single-valley 
semiconductors,  has  been  extended  to  multi-valley  semiconductors.  A  two-valley 
model  has  been  implemented  and  the  results  compared  to  those  obtained  using 
Monte  Carlo  methods.  Excellent  agreement  has  been  obtained  in  situations  tiriieie 
a  single  valley  model  is  inadequate  due  to  significant  inter-valley  carrier  transfer.  The 
macro-kinetic  model  is  considerably  easier  to  implement  numerical^  than  those  based 
on  the  Boltzmaim  transport  equation  and  requires  orders  of  magnitude  less  computer 
time  to  nm.  The  theory  is  presented  in  Sec.  IIL 

(2)  The  theory  described  above  has  been  successfully  applied  to  modeling  of 
electron  transport  in  GaAs  subjected  to  rapidly  varying  fields.  This  scheme  provides 
an  efficient  scheme  for  device  modeling,  and  overcomes  the  difficulty  of  evaluating 
the  relaxation  times  in  multi-valley  conservation  equations  without  a  priori 
assumption  of  a  displaced-Maxwellian  distribution  for  the  carriers.  This  application 
is  discussed  in  Sec.  IV. 

(3)  The  modelling  of  charge  transport  in  GaAs  requires  a  number  of  material 
parameters,  such  as  the  effective  mass,  acoustic  deformation  potential,  various 
phonon  energies,  energy  separation  between  valleys,  inter-valley  and  intra-valley 
deformation  potentials,  etc.  At  present,  there  are  large  imcertainties  in  the 


experimental  and  theoretical  values  for  some  of  these  parameters.  The  Monte  Carlo 
code  developed  in  the  previous  funding  period  has  been  modified  to  develop  a  set 
of  material  parameters  for  GaAs  by  minimizing  the  error  between  calculated  and 
measured  values  of  the  transport  parameter  over  the  range  of  applied  electric  fields 
for  which  experimental  results  are  available.  This  investigation  and  the  results 
obtained  are  presented  in  Sec.  V. 

(4)  The  experimental  fiicility  to  investigate  the  dynamics  of  field  collapse  in  high 
power  photoconductive  switches  has  been  completed  and  preliminary  results  have 
been  obtained.  A  disciission  of  the  facility  is  presented  in  Appendix  A. 
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A  number  of  presentations  were  also  made  during  this  period.  Abstracts  are 
given  in  Appendix  B.  Also  during  this  period,  one  student  (M.  Cheng)  received  his 
Ph.D.  and  is  currently  on  the  faculty  at  the  University  of  New  Orleans.  A 


postdoctoral  fellow  (A.  Garvie)  joined  the  research  and  has  been  responsible  for 


setting  up  the  experimental  facih'ty. 
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DDL  A  Theaiy  of  Nonequilibriiim  Cankr  lYansport  in  Multivalley  Semiconducton 


I.  INTRODUCTION 

The  observation  of  the  Gunn  effect confirmed  the 
prediction^  *  of  electron  transfer  between  light-mass  and 
heavy-mass  valleys  in  GaAs  and  InP.  This  transfer  has  sub¬ 
sequently  been  found’  *  to  have  a  significant  influence  on 
electron  transport  in  multivalley  semiconductors,  particu¬ 
larly  in  nonequilibrium  conditions  and  also  at  fields  where 
intervalley  scattering  is  a  dominant  process.  For  example,  in 
high-frequency  devices  with  a  short  active  channel,  the  in¬ 
tervalley  transitions  play  an  important  role  near  the  cathode 
contact.’  This  can  generally  be  said  of  any  device  with  small 
temporal  and/or  spatial  scale.  Similarly,  the  velocity-under¬ 
shoot  behavior*  of  electrons  in  GaAs  subjected  to  a  picose¬ 
cond  pulse  electric  field  is  explained,  in  part,  in  terms  of  the 
electron  transfer  from  upper  valleys  to  the  central  valley, 
which  scattering  process  results  in  a  randomization  of  the 
momentum. 

At  the  kinetic  level,  the  semiclassicai  description  of  car¬ 
rier  dynamics  in  multivalley  semiconductors  is  given  in 
terms  of  the  phase-space  distribution  functions  for  the  carri¬ 
ers  in  each  valley.  These  distributions  obey  Boltzmann  trans¬ 
port  equations  (BTEs)  that  arc  coupled  by  intervalley  scat¬ 
tering  processes.  Difficulties  in  dealing  with  the  scattering 
terms  in  the  coupled  system  of  equations  make  this  approach 
unattractive.  To  simplify  the  problem,  these  terms  have  been 
approximated  by  using  the  relaxation  approximation  for 
both  intervalley  and  intravalley  scatterings^  * ;  however,  this 
approach  loses  its  generality  at  high  fields.  An  alternate  pro¬ 
cedure  to  the  BTE  for  obtaining  the  distribution  function  in 
each  valley  is  to  use  Monte  Carlo  (MC)  methods.  Although 
easy  to  implement,  this  procedure  is  computationally  time 
consuming,  particularly  for  device  simulations. 

A  number  of  approaches  have  been  proposed  with  the 
objective  of  obtaining  a  model  that  is  easier  to  implement 
while  including  the  pertinent  physics.  Most  of  them  make 
use’"’*  of  a  finite  set  of  moments  of  the  BTE.  The  difficulty 
in  implementing  this  approach  is  the  determination  of  the 


unknown  coefficients  appearing  in  the  truncated  moment 
equations.  This  is  especially  true  in  multivalley  semiconduc¬ 
tors,  which  require  a  set  of  moment  equations  for  each  val¬ 
ley.  The  difficulty  has  resulted  in  the  fact  that  most  treat¬ 
ments  of  transport  in  multivalley  semiconductors  using 
moment  equations  assume  only  a  single  valley.  In  this  ap¬ 
proximation,  a  single  distribution  function,  averaged  over 
the  valleys,  is  taken  to  represent  the  distributions  in  all  val¬ 
leys;  consequently,  the  moments  correspond  to  averages 
over  the  vaUeys.  Moreover,  with  this  distribution,  the  un¬ 
known  coefficients  in  the  “equivalent”  single-valley  moment 
equations  can  be  easily  evaluated.*  However,  the  dynamics 
of  carrier  transfer  arc  not  properly  taken  into  account  when 
ining  this  single-valley  description.*'”''*  As  an  alternative, 
and  in  order  to  keep  a  multivalley  description,  the  assump¬ 
tion  is  made  that  the  distribution  function  in  each  valley  is  a 
(different)  dbplaced  Maxwellian  distribution.’  '’''*  These 
distributions  are  then  used  to  evaluate  the  unknown  coeffi¬ 
cients  in  the  moment  equations  for  each  valley.  However, 
this  assumption  is  rather  severe,  since  the  distribution  may 
be  appreciably  different  from  a  displaced  Maxwellian  in 
some  semiconductors. 

In  this  paper  a  theory  of  nonequilibrium  carrier  trans¬ 
port  in  multivalley  semiconductors  is  presented.  It  is  based 
on  the  concept  of  a  macrokinetic  distribution  ( MKD ) 
and  on  moments  of  the  BTE.  No  a  priori  assumptions  are 
made  in  this  theory  about  the  shape  of  the  carrier  distribu¬ 
tion  in  momentum,  k,  space.  The  MKD  is  a  momentum- 
dependent  function  with  macroscopic  space-time  scales  and 
obeys  an  equation  that  is  macroscopically  equivalent  to  the 
BTE.  The  MKD  can  also  be  obtained  using  MC  methods. 
The  unknown  parameters  and  rates  appearing  in  the  multi¬ 
valley  moment  equations  can  be  evaluated  using  the  MKD. 
It  is  shown  that  in  the  case  of  strong  intervalley  transitions, 
the  single-valley  description  does  not  account  for  energy  and 
momentum  exchanges  between  valleys.  This  leads  to  a  sig¬ 
nificant  difference  between  the  multi-  and  single-valley  be¬ 
havior  for  the  relaxation  of  mean  energy  and  momentum.  '* 
For  simplicity,  a  semiconductor  test  model  is  set  up  to  devel¬ 
op  the  theoretical  details  of  the  macrokinetic  model  and  to 


carry  out  the  numerical  calculations.  The  band  structure 
and  the  scattering  mechanisms  used  in  this  model  are  given 
in  the  Appendix.  Results  are  presented  for  the  response  of 
electrons  to  a  step  change  in  electric  field  and  compared  to 
those  obtained  using  MC  methods  and  an  equivalent  single¬ 
valley  model. 

In  the  next  section,  a  •'rief  overview  of  the  macrokinetic 
theory  for  a  single-valley  c,and  is  given.  The  multivalley  de¬ 
scription  is  shown  in  Sec.  III.  In  Sec.  IV  the  results  and 
comparisons  are  presented.  Concluding  remarks  are  given  in 
Sec.  V. 

II.  THE  MACROKINETIC  MODEL 

This  model  has  been  discussed  for  a  single-valley  band 
in  a  previous  paper'* ;  consequently,  only  a  brief  description 
is  presented  here.  It  is  a  prior/  assumed  that  for  the  situation 
of  interest  a  finite  set  of  moment  equations  can  be  used  to 
describe  carrier  behavior.  A  MKD,/^# ,  is  introduced  that  is 
macroscopically  equivalent  to  the  microscopic  phase-space 
distribution, /(k,r,r).  The  obeys  an  equation'*  that  is 
equivalent  to  the  BTE  in  the  space-time  scale  of  the  moment 
equations.  The  equation  for  and  the  moment  equations 
constitute  the  model.  This  set  can  be  used  to  describe  the 
nonequilibrium  behavior  of  carriers  in  the  time  scale  corre¬ 
sponding  to  the  characteristic  times  in  the  moment  equa¬ 
tions. 

Since  r  >  r,  >  in  semiconductors,  where  r,  r,,  and 
T„  are  carrier  density,  energy,  and  momentum  relaxation 
times,  respectively,  the  moment  equations  can  be  ordered 
and  a  finite  set  selected  depending  on  the  level  of  description 
or  degree  of  space-time  coarseness  desired  in  the  model.  To  a 
given  finite  set  of  moment  equations  corresponds  a  MKD 
that  serves  to  close  the  set.  The  equation  for  this  MKD  is 
obtained  from  the  BTE.  For  a  single-valley  band,  the  BTE  is 
written  as 

dj  +  +  (1) 

n  I  otic 

where  v  is  the  velocity  and  E  is  the  electric  field.  The  corre¬ 
sponding  first  three  moment  equations  for  the  single  valley 
are 

d,n  -I-  V*(nu)  =  n/r,  (2a) 

(ne) -t- V*<eT>  —  ^nE-u  =  —  [n(?  -  tg  )/r,  ],  (2b) 

^,(np) -I- V<p-v) -gnE=  —  (np/r„),  (2c) 

where  the  bracket  represents  an  average  over  the  distribu¬ 
tion,  n  is  the  carrier  density,  u  is  the  average  velocity,  ?  is  the 
mean  energy, p  is  the  average  momentum,  and  fg  >s  the  ther¬ 
mal  energy.  The  carrier  density,  energy,  and  momentum  re- 


taxation  times  are  defined  by 

7  = 

(3a) 

it(f-eo) 

(3b) 

(3c) 

Subsequently,  the  various  levels  of  description  that  can  be 
obtained  using  Eqs.  (2)  are  discussed. 


A.  The  Si  description,  valid  for  time  scale 

In  this  time  scale,  the  dynamics  of  the  carriers  can  be 
described  by  Eq.  (2a).  Consequently,  the  corresponding 
MKD  is  taken  to  have  the  same  space-time  dependence  as 
the  carrier  density;  that  is, 

/(k,r,r)  =/JJ>[k,n(r,r)].  (4) 

Using  Eqs.  (2a)  and  (4)  in  Eq.  ( 1 ),  the  equation  for is 
obtained.  Solving  for Z^’  and  using  it  to  evaluate  u  and  r  in 
Eq.  (2a)  leads  to  the  well-known  drift-diffusion  equation  for 
n(r,r).  That  is, 

nu  =  npE  — DVn,  (5) 

where  p  is  the  mobility  and  D  is  the  diffusion  coefficient  and 
are  obtained  in  terms  of  the  known  function ZU’-'^ 

B.  The  Sa  description,  valid  for  time  scale  ^r. 

In  this  case,  the  description  is  in  terms  of  a  MKD  that 
depends  on  the  carrier  density  and  mean  energy, 

Z[k,r,r)  =ZiJ’[k,n(r,r),e(r,r)],  (6) 

and  the  corresponding  first  two  moment  equations  [Eqs. 
(2a)  and  (2b)  ].  The  equation  for Z«’  has  been  shown'*  to 
be,  for  uniform  conditions, 

=  (7a) 

where 

qE^  =  [?(r)  -  fg  ]/u,r„  (7b) 

with 

nu,  =  j/ji'(k,n,e)7dk.  (7c) 

Equation  (7a)  represents  a  system  that  obeys  a  steady-state 
BTE  with  an  equivalent  field,  E^ ,  as  the  driving  force,  u,  is 
the  “steady-state”  average  velocity  that  is  seen  to  depend  on 
the  time- varying  mean  energy.  Equation  (7c)  implies  that  in 
this  description  the  mean  momentum  of  carriers  responds 
instantaneously  to  an  equivalent  field,  . 

The  unknown  rates  in  Eqs.  (2a)  and  (2b)  can  be  ob¬ 
tained  in  terms  of  the  mean  energy  by  evaluating  the  inte¬ 
grals  in  Eqs.  (3a)  and  (3b)  using  the  MKD  obtained  from 
Eq.  (7a).  The  Sj  MKD  can  be  obtained,  for  example,  from 
MC  simulations  and  the  rates  tabulated  as  functions  of  the 
mean  energy  or  equivalent  field.  For  a  given  problem,  it  is 
then  only  necessary  to  solve  Eqs.  (2a)  and  (2b). 


C.  The  Sj  description,  valid  for  time  scaie  ^  r„ 

The  system  of  equations  that  describe  the  system  in  this 
time  scale  consists  of  the  lint  three  mooMBt  eqoatioiis  [  Eqs. 
(2a)-(2c)]  oftheBTEandthecorrespoodiag  MKD,  which 


is  taken  to  depend  on  carrier  density,  mean  energy,  and  aver¬ 
age  momentum: 

/(k,r,f)  =/M’[k,n(r,r),?(r,r),k(r,r)]. 

In  a  previous  paper, a  relaxation  process  is  used  to  model 
the  evolution  of  the  smaller  time-scale  distribution  (/«’) 
into  the  larger  time-scale  distribution  ( /j^’ ) .  An  alternative 
is  to  approximate with in  the  evaluation  of  the  aver¬ 
aged  rates  [Eqs.  (3a)-(3c)  ].  This  simple  approach  is  valid 
when  the  momentum  relaxation  time  is  relatively  small  com¬ 
pared  to  the  energy  relaxation  time.  It  is,  in  general,  true  for 
GaAs  and  Si  unless  the  situation  is  considered  at  relatively 
low  fields. Similarly  to  the  S,  and  descriptions,  the 
unknown  coefficients  in  Eqs.  (2a)-(2c)  can  be  evaluated 
using/iJ>. 

III.  THE  MULTIVALLEY  DESCRIPTION 

For  simplicity  a  two-valley  model  is  used  to  illustrate 
the  theory.  The  BTEs  for  this  model,  assuming  spatial  uni¬ 
formity,  are  given  by 


where /  (kj,  is  the  distribution  function  in  the  rth  valley, 

-J/>„(k„k;)y;(k,)dk;, 

® 

\  dt  /  {.inter  J 

(fL 

with  ij  —  1,2.  The  transition  rates  per  state,  are  given  in 
the  Appendix.  The  first  term  on  the  right-hand  side  of  Eqs. 
(8)  represents  the  electron-phonon  collisions  inside  the 
same  valley,  and  the  last  two  terms  represent  the  intervalley 
transitions. 

The  conservation  { macroscopic )  equation  correspond¬ 
ing  to  a  microscopic  physical  quantity,  d,  (k, ),  defined  in 
each  valley  can  be  obtained  by  averaging  over  the  distribu¬ 
tion  and  using  the  BTE.  The  resulting  equation  is 


=  (£^)  , 

'  dt  /  inlri  V  dt  /  {.inter 

where 

\  dt  J  /.inter  \  dt  J  /.inter 

■ 

\  dt  /  /.inter 

and  {(^)  =  ri^  —  n,^i  +  Pof  =  1*  ^i»  “d  fit® 
moment  equations  for  density,  average  momentum,  and 
mean  energy  are  obtained. 

A.  Density  conservation  equations, 

<9,n,  =/i2V2, -rt,v,j,  (9a) 

=n,v,2 -/ijVji,  (9b) 

where 

=  J±l) 

\  dt  /  l.inio'  V  dt  Junta 

=  Jp,2(k,,k2)/,  (k,  )rfk,  rfk,,  (10a) 

'  '  V  *  V  *  Junta 

=  J  Pj.  (k2.k,  )/2  (kj  )</k,  rfk2.  ( 10b) 

Therefore,  the  total  intervaUey  transition  for  carrier  density 
in  the  multivalley  system  vanishes.  That  is, 

fil)  =(i!L]  =0, 

\dt/  inter  \dtJt  .inter  \dtJ  24ater 

where  n  =  n,  -f  ^2.  This  implies  that  in  the  S,  regime, 
where  the  time  scale  of  interest  is  much  greater  than  energy 
and  momentum  relaxation  times,  the  single-valley  descrip¬ 
tion  is  sufficient  to  determine  carrier  behavior  in  multivalley 
semiconductors. 


B.  Momentum  conservation  equation,  ^,(k,)=;fik,  =;p, 
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with 
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J  P^(k,,k^)dk^, 
J  />„(k,.k;)(/k;. 


T“ '  and  T~  given  in  the  Appendix,  are  the  intra valley  and 
intervalley  scattering  rates,  respectively.  r„„  and  r„,j  are  the 
intravalley  and  intervalley  momentum  exchange  times,  re¬ 
spectively.  The  total  intervalley  momentum  exchange  rate  in 
this  model  is  given  by 


where  <p>  =  np  =  /i,p,  n^2^^Pi  =  In  the  case  of 

strong  intervalley  transitions,  t,J  '  is  comparable  to  and 
thus  from  Eqs.  ( 12a)  and  ( 12b)  r ‘J  is  comparable  to  t“„‘. 
Consequently,  the  total  intervalley  momentum  exchange 
rate  [Eq.  ( 13)  ]  is  comparable  to  the  total  intravalley  mo¬ 
mentum  exchange  per  unit  time  which  is  given  by 


^  "iPi  "jPi 
\  dt  Jmtn  r„21 

Note  that  the  single-valley  approximation  does  not  properly 
account  for  the  momentum  exchange  between  valleys. 


C.  Energy  conservation  equation,  ^,(k;)=c 
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where  the  intravalley  and  intervalley  energy  exchange  times, 
r,„  and  respectively,  are  defined  as 


n,(€,  -  Co)  _ 

V  /  infra 

=  ft/ k'f/ik, 

-  [“k,  . 
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The  phonon  energy,  Aty^, ,  and  the  energy  separation  between 
valleys,  are  given  in  the  Appendix.  It  is  noted  that 


’■<121  7^ ’■<221  and  7-£,,2  ^^’■<212-  ^ven  in  the  case  where  Ae^, 
equals  zero  since  ^0.  The  total  energy  exchange  per 
unit  time  between  valleys  is 

(d{€)\  n2f2  «2f2  ,  n,?,  _ 

I— — 1  - - -I- -  (16a) 

V  Ol  /inter  ■’■<121  ’’<221  ’’<212  ’’<112 
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=^^21  («2V2I  -«lV,2) 

\  Ot  /inter 

-l-Aa),j/i,(v,2„4 -v,j,„),  (16b) 

where  )  is  the  averaged  transition  rate  from  the  fth 

to /th  valley  through  the  process  of  phonon  emission  (ab¬ 
sorption),  and  given  by  Eq.(  10)  with  only  the  phonon  emis¬ 
sion  ( absorption )  transition  rate.  The  first  term  on  the  right- 
hand  side  RHS  of  Eq.  ( 16b)  equals  Ae2i  <^r'>i  >  and  does  not 
vanish  for  Afji  ^0  unless  carriers  are  in  equilibrium  be¬ 
tween  valleys.  This  clearly  shows  that  for  semiconductors 
with  Af,y  7^0  the  mean  energy  can  only  reach  steady  state  in 
the  same  scale  as  the  carrier  density  in  each  valley.  The  sec¬ 
ond  and  third  terms  on  the  RHS  of  Eq.(  16b)  are  significant 
in  comparison  to  the  intravalley  e&ergy  exchange  when  the 
intervalley  transition  is  dominant  Again,  the  dynamics  of 
energy  exchange  between  valleys  is  not  pr(^)erly  taken  into 
account  when  using  a  single-valley  description. 

To  obtain  the  carrier  density,  mean  energy,  and  momen¬ 
tum  from  the  multivalley  moment  equations,  the  (un¬ 
known)  averaged  rates  [Eqs.  (9),  ( 1 1 ),  and  ( 14)  ]  have  to 
be  determined.  To  do  so  and  proceeding  as  for  a  single- valley 
band  (Sec.  II),  the  MKD  in  the  rth  valley, /m,,  is  used  to 
directly  evaluate  Eqs.  ( 10),  ( 12),  and  (IS)  (this  is  not  nec¬ 
essary  in  practice  for  a  single-valley  band  since  the  desired 
rates  can  also  be  calculated  using  the  steady-state  equations 
if  the  steady-state  moments  are  known’’ ).  Depending  on 
the  time  scale  of  interest,  the  MKD  and  the  corresponding 
macrokinetic  equation  in  each  vaUey  can  be  obtained  for 
each  level  of  description  as  shown  in  Sec.  II.  In  the  next 
section,  results  are  presented  using  an  S3  description  (in  the 
momentum  scale)  with  an  Sj  MKD,  /]^/,  i.e.,  with 
/  Mi  ~fMi  •/«/  can  be  obtained  in  terms  of  mean  energy  from 
MC  methods  and  the  integrals  [Eqs.  (10),  (12),  and  (IS)] 
can  be  similarly  determined.  The  resulting  (closed)  set  of 
moment  equations  constitute  a  mode!  of  nonequilibrium  car¬ 
rier  transport  in  multivalley  semiconductors. 

(V.  RESULTS 

Results  are  presented  in  this  section  for  the  band  model 
discussed  in  Sec.  III.  This  band  model  has  two  different  val¬ 
leys  with  effective  masses  equal  to  O.OSmg  (lower  valley) 
and  0. 3mo  ( upper  valley ) ,  respectively,  where  mo  is  the  free- 
electron  mass.  The  e-k  relation  of  carriers  in  each  valley  is 
characterized  by 

e,  =  A^JIc  ,V2m, . 

The  results  obtained  with  this  model  are  also  compared  to 
those  obtained  using  a  single-valley  “equivalent”  descrip¬ 
tion'^  and  multi  valley  MC  simulations.'^  For  further  sim- 


plilication,  a  test  model  is  also  used  for  the  k-space  transition 
rates.  The  scattering  model  is  given  in  the  Appendix.  These 
transition  rates  are  used  in  all  the  calculations  presented  in 
this  paper. 

To  illustrate  the  dynamics  of  the  carriers  in  multivalley 
semiconductors,  the  transient  behavior  of  electrons  (in  the 
test  model )  influenced  by  a  step  change  in  electric  field  is 
presented.  The  applied  field  is  along  the  z  direction  and  equal 
to  £■ ,  for  /  <  r  I  and  £2  for  r>r , .  For  /  <  / , ,  the  system  is  in 
steady  state.  A  finite  difference  technique^^  has  been  used  to 
solve  the  moment  equations  for  the  single-valley  descrip¬ 
tions  [Eqs.  (2a)-(2c)  without  space  dependence]  and  for 
the  multivalley  description  [Eqs.  (9),  (11),  and  ( 14)  ]. 

In  the  first  example,  £,  =  5  kV/cm,  £2  =  15  kV/cm, 
and  r  I  =  0. 3  ps.  The  behavior  in  time  of  the  total  mean  ener¬ 
gy  and  average  velocity  obtained  from  numerical  solutions 
in  Eqs.  (9),  ( 1 1 ),  and  ( 14)  are  shown  in  Fig.  1.  The  rates  in 
these  equations  are  evaluated  from  the  integrals,  Eqs.  ( 10), 
( 12),  and  (IS)  using  the  MKD  in  each  valley.  The  MKD  is 
obtained  from  MC  simulations  in  term  of  energy,  as  dis- 


FIG.  1  Evolution  of  total  mean  energy  and  average  velocity  as  a  conse¬ 
quence  of  a  step  change  in  applied  field  from  5  to  15  kV/cm.  (a)  Mean 
energy  and  (b)  average  velocity. 


cussed  in  Sec.  II.  The  results  obtained  using  other  methods 
are  also  shown.  It  can  be  seen  that  the  single-valley  variables 
respond  to  the  change  in  electric  fields  faster  than  those  cor¬ 
responding  to  a  multi  valley  model  (both  the  moment  de¬ 
scription  and  the  MC  method).  Moreover,  the  transient  en¬ 
ergy  and  velocity  obtained  from  the  multivalley  moment 
equations  based  on  the  macrokinetic  model  are  in  excellent 


FIG  2.  Evolution  of  normalized  populatian.  mean  energy,  and  average  ve¬ 
locity  of  the  camera  in  each  valley  as  a  conaeqneiKe  of  a  step  change  in 
applied  field.  Same  conditions  as  Fig.  I.  (a)  Population,  (b)  mean  energy, 
and  (c)  average  velocity. 


FIG.  3.  Evolution  of  total  mean  energy  and  average  velocity  of  the  test  mod¬ 
el  as  a  consequence  of  a  step  change  in  applied  field  from  1 3  to  3  k  V  /cm.  ( a ) 
Mean  energy  and  (b)  average  velocity. 


agreement  with  those  obtained  from  MC  calculations  (see 
Fig.  1). 

The  slower  carrier  relaxation  observed  using  a  multival¬ 
ley  description  as  compared  to  a  single-valley  description 
can  be  explained  with  the  aid  of  Eqs.  (13)  and  (16b).  The 
slowest  time  scale  in  the  model  is  that  of  the  carrier  density. 
The  slow  response  of  the  carrier  density  to  changes  in  the 
field  (see  Fig.  2(a)]  influences  the  change  in  mean  energy 
through  the  first  term  on  the  RHS  of  Eq.  (16b),  and  thus  the 
velocity  response  since  the  moment  relaxation  times  are 
strongly  dependent  on  energy.  In  the  case  of  strong  interval¬ 
ley  transitions,  the  intervalley  momentum  and  energy  ex¬ 
changes  in  Eq.  (13)  and  the  last  two  terms  in  Eq.  ( 16b)  also 
determine  the  dynamical  evolution  of  energy  and  momen¬ 
tum.  The  evolution  of  carrier  density,  mean  energy,  and 
average  velocity  in  each  valley  is  shown  in  Fig.  2.  This  also 
shows  excellent  agreement  between  the  macrokinetic  model 
using  a  multivalley  description  and  the  MC  method. 

Similar  comparisons  are  shown  in  Figs.  3  and  4  for  a  step 
change  in  field  from  E,  =  15  kV/cm  to  Ej  =5  kV/cm  and 


Time  (pa«c) 


FIG.  4.  Evolution  of  normalized  population,  mean  energy,  and  average  ve¬ 
locity  of  carriers  in  each  valley  as  a  consequence  of  a  step  change  in  applied 
field.  Same  conditions  as  Fig.  3.  (a)  Population,  (b)  mean  energy,  and  (c) 
average  velocity. 


f ,  =  0.5  ps.  The  slower  responses  of  mean  energy  and  aver¬ 
age  velocity  to  the  change  in  field  using  multivalley  descrip¬ 
tion  as  compared  to  a  single-valley  description  (Fig.  3)  can 
also  be  seen.  This  is  due  to  the  slower  carrier  density  re¬ 
sponse  [see  Fig.  4(a)  ]  to  the  time-varying  field  as  discussed 
for  the  step-up  change.  The  evolution  of  carrier  density, 
mean  energy,  and  average  velocity  in  each  valley  using  the 
macrokinetic  model  and  the  MC  method  is  shown  in  Fig.  4. 
Again,  excellent  agreement  is  obtained  between  these  two 
descriptions. 

For  each  step  field,  MC  simulations  ( 50  000  electrons 
have  been  used)  require  more  than  4  h  of  computation  time 
in  the  Micro-Vax  II.  However,  solving  the  moment  equa¬ 
tions  for  the  single- valley  description  and  for  the  multivalley 
description  takes  about  1.7  and  7.0  s,  respectively,  in  the 
same  computer. 

V.  CONCLUDING  REMARKS 

A  theory  of  nonequilibrium  transport  of  carriers  in  mul¬ 
tivalley  semiconductors  has  been  presented.  It  is  based  on 
the  macrokinetic  model for  the  distribution  function  in 
each  valley.  It  is  shown  that  in  the  case  of  nonequUibrium 
and  strong  intervalley  scatterings  the  transfer  effect  of  carri¬ 
ers  between  valleys  is  not  properly  taken  into  account  when 
using  a  single-valley  description.  The  carriers  are  observed 
to  respond  faster  to  changes  in  the  applied  field  when  using  a 
single-valley  “equivalent"  model  as  compared  to  a  multival¬ 
ley  model.  That  is,  the  single-valley  model  by  not  properly 
accounting  for  the  dynamics  of  intervalley  energy  and  the 
momentum  exchanges  drives  the  carrier  to  equilibrium  fas¬ 
ter  than  the  multivalley  model.  Moreover,  since  the  scales  at 
which  the  relaxation  occurs  is  that  of  the  moment  equations, 
the  MKD  can  correctly  yield  the  macroscopic  scattering 
rates.  This  can  be  seen  from  the  excellent  agreement  between 
the  results  from  multivalley  moment  model  and  the  MC  re¬ 
sults. 
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APPENDIX 

The  multivalley  semiconductor  test  model  used  in  this 
paper  has  a  central  valley  with  mass  equal  to  O.OSmg  and  two 
equivalent  satellite  valleys  with  mass  equal  to  O.OJmQ.  The 
energy  difference  between  the  lower  ( central )  valley  and  up¬ 
per  valleys  is  Af,,  =0.15  eV.  All  valleys  have  a  parabolic 
band  structure.  The  e-k  relation  is  given  by  e  =  ^k  ]/2m, . 

The  transition  rate  for  intravalley  scattering  from  state 
k,  to  k;  is  given  by 

S„(k,,k,')  =  (2ff’/fi)|//,., |’<5(e,'  -  e,  ±  Au,)  (Al) 

and  by 

5„(k,,k,)  =  (2ff-/^)!//„;’^(f,  -  f,  ±fuo„  -  Af,,)  (A2) 


for  intervalley  scattering  from  state  k,  to  k,,  where 
ij  =  1,2,/^y,  Af|2  =  —  Ae,| ,  and  H,j  (or  //,, )  is  the  ma¬ 
trix  element  between  the  initial  and  the  final  states.  Thus, 
transition  rate  per  state  in  Sec.  II  is  given  by 


P„(k,,k;)  =  (F/27r)5„(k„k:), 

/>,^(k„k,)  =  (F/2;r)5,(k„k,), 

where  V  is  the  crystal  volume.  The  scattering  mechanism  in 
each  valley  is  given  as  follows; 

(i)  Intravalley  scattering;  The  intra valley  phonon  ener¬ 
gy  in  each  valley  is  taken  to  be  equal,  and 
Au,  =fuo2  =0.035  eV.  The  matrix  element,  is  as¬ 
sumed  to  be  constant  for  either  the  emission  process  or  the 
absorption  process.  Therefore,  Eq.  ( 1 )  can  be  rewritten  as 


5„(k„k;)  5, 

n 


(A3) 


iv,5(e;  -  e,  -  fko, ) 

(N,  +  l)5(e;  -€,+fio},)\ 

where  is  a  constant  and  N,  is  the  phonon  distribution.  In 
this  model,  N,  is  assumed  to  be  in  equilibrium  and  given  by 


exp(fko,/kgT)  —  1 

where  kgT  is  the  thermal  energy.  In  Eq.  (A3)  A,  and 
A,  +  1  correspond  to  the  processes  of  phonon  absorption 
and  emission,  respectively.  The  total  scattering  rate  out  of 
the  state  k^  for  this  process  is,  therefore,  given  by 


1 

r,v(*,) 


5„(k„k:)dk;. 


For  a  parabolic  valley,  the  total  scattering  is  found  to  be 


(A4) 


where 


'  itH* 

e;  =  e,  ±  Au,. 

In  this  test  model,  D,  =2.5x10'^  s  '  eV " and 
Dj  =5.0X10’’  s  '  eV-'^^ 

(ii)  Intervalley  scattering;  Intervalley  phonon  energy  is 
taken  to  be  Ai;,2  =  Auji  =  0.035  eV.  Similarly,  Hj,  is  taken 
to  be  constant.  Equation  ( A2)  is  rewritten  as 


1  ( A,^  +  1  )fi(fy  -  e,  +  ^,J  -  ) 


.  (A5) 


where  B,j  is  constant,  A,y  is  the  intervalley  phonon  distribu¬ 
tion  which  is  in  equilibrium  and  given  by 


''  exp(fko,/kgT)  —  \ 

In  Eq.  ( A5)  A,y  and  A^  -|-  1  correspond  to  the  processes  of 
phonon  absorption  and  emission,  respectively.  Similarly,  the 
total  scattering  rate  out  of  k,  for  this  process  is  written  as 
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Tyfk,) 


(2fl-) 


-  J  \(k„k^)</ky. 


For  our  model,  it  can  be  rewritten  as 


IV.  A  Transport  Model  for  Carriers  in  Small  SpatiaVTemporal  Scale  Cnmjyninri 
Semiconductor  Devices 


••t 

I.  Introduction 

>  '  Device  modeling  requires  a  proper  transport  model  to  describe  the  carrier  transport 
phenomena  that  characterize  device  performances,  such  as  the  upper  frequency  limit  and 
power  efficiency.  The  choices  of  transport  models,  in  general,  depend  on  space/time  scales 
'of  interest.  For  example,  the  conventional  model  for  device  analysis  is  based  on  th^  drift- 
diffusion  model  which  assumes  a  quasi-equilibrium  distribution  of  carriers  and  thus  can  not 
account  for  non-equilibrium  behavior  due  to  rapidly  space/time-varying  fields,  such  as 
overshoot  [1-5]  and  undershoot  [6]  phenomena.  This  model  is,  therefore,  inadequate  when 
frequencies  approach  several  10  GHz  and/or  devices  are  scaled  down  to  submicrom^er  size. 
In  such  small-scale  devices,  a  transport  model  including  pertinent  physics  of  non-equilibrium 
phenomena  is  thus  desirable.  Most  of  transport  models  that  have  been  used  [7-13]  to  account 
for  non-equilibrium  phenomena  of  carriers  in  compound  semiconductor  devices  is  based  on 
the  single-valley  conservation  equations.  It  has  been  assumed  that  the  energy  and  momentum 
relaxation  times  in  the  conservation  equations  only  depend  on  the  carrier  mean  energy 
(theoretically,  on  all  moments  of  the  carrier  distribution  fonction  [14,15])  and  are  described 
by  the  phenomenological  equations  [7,8].  In  compound  semiconductor  devices,  carrier 
behavior  is  described  by  the  distribution  function  in  each  valley.  These  distributions  in 
different  valleys  obey  the  Boltzmann  transport  equations  (BTE’s)  that  ate  coupled  by  inter- 
valley  scattering  processes.  Consequently,  the  macroscopic  description  is  given  by  a  set  of 
multi-valley  conservation  equations  which  are  thus  coupled  by  carrier  inter-valley  transitions. 
However,  the  relaxation  times  in  the  multi-valley  conservation  equations  can  not  be 
determined  from  the  phenomenological  equations  [7,8].  The  difficulty  in  determining  these 
times  has  led  to  the  fact  that  the  use  of  the  conservation  equations  for  carrier  dynamics  in 
compound  semiconductor  devices  has  assumed  [1,2,7-13]  only  a  single  valley  where  the 


carrier  density,  mean  energy  and  average  momentum  are  averaged  over  all  valleys. 

Inter-valley  transfer  has  significant  influences  on  non-equilibrium  carrier  dynamics  in 
compound  semiconductor  devices.  It  has  been  previously  found  [14-16]  that  the  use  of  the 
single-valley  conservation  equations  for  the  carrier  dynamics  in  GaAs  results  in  faster 
response  to  a  rapid  change  in  electrical  field  than  actually  occurs.  This  behavior  is  caused 
by  the  slow  averaged  inter-valley  transitions  [16]  that  are  not  properly  taken  into  account  in 
the  single-valley  conservation  equations.  In  situations  where  inter-valley  scattering  is 
dominant,  the  multi-valley  conservation  equations  are  therefore  necessary.  In  order  to 
determine  the  relaxation  times  for  solving  the  multi-valley  conservation  equations  and  to 
account  for  inter-valley  transfer  effects,  the  assumption  has  been  made  in  most  papers  [17- 
23]  that  the  distribution  function  in  each  valley  is  a  displaced  Maxwellian.  These 
distributions  are  then  used  to  evaluate  the  relaxation  times.  However,  it  is  well-known  that 
in  hot  electron  situations  the  electron  distribution  function  appreciably  differs  from  a 
displaced  Maxwellian  in  compound  semiconductors  [24,25]. 

In  this  paper,  a  non-equilibrium  transport  model  for  compound  semiconductor  devices 
has  been  developed  based  on  the  concept  of  the  Macro-Kinetic  Distribution  (MKD)  functions 
[14-16]  (whose  time  scales  are  characterized  by  moments  of  the  distribution).  The  model 
provides  an  efficient  scheme  for  device  modeling,  and  can  account  for  the  transfer  effects 
without  a  priori  assumption  of  a  displaced  Maxwellian  distribution.  The  mathematical 
formulation  of  this  multi-valley  transport  model  is  given  in  the  next  section.  The  resulting 
model  is  applied  to  electron  transport  in  three-valley  GaAs  with  non-parabolic  bandstructure 
subjected  to  rapidly  time-varying  fields.  The  results  compared  to  those  obtained  from  a 
Monte  Carlo  meth^  [25]  are  presented  in  Sec.HI.  Conclusions  are  given  in  Sec.IV. 

1 

n.  The  Multi-Vallev  Transport  Model 


The  exact  description  of  carrier  transport  in  compound  semiconductors  subjected  to 
space/time-varying  fields  in  semi-classical  regime  is  given  by  the  BTE’s  that  are  coupled  by 
inter-valley  scattering  processes  and  written  as 


.il. 


V,= 


dt) 


inm 


(1) 


ifiHUr 


where /jYriATj  ,t)  is  the  ith-valley  distribution  in  phase  space,  is  the  number  of  valleys,  w,- 
is  the  carrier  velocity,  k,  is  the  wave  vector,  E  is  the  electric  field,  and  ft  is  the  Planck 
constant  divided  by  2t.  TTie  first  term  of  the  RHS  represents  the  mtra-valley  scattering,  and 
the  inter-valley  scattering  processes  out  of  and  into  the  ith  valley  are  describe  by  the  second 
and  third  terms  of  the  RHS,  respectively.  The  carrier  density,  energy  and  n'  amentum 
conservation  equations  in  each  valley  can  be  obtained  from  Eq.  (1)  by  averaging  over  the 
distribution,  and  are  given  by 
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where  kg  is  the  Boltzmann  constant,  ind  is  the  electron  temperature  related  to  carrier 
mean  energy  by 

m*  is  the  macroscopic  effective  mass  of  carriers,  and  n, ,  ^ ,  and  U;  are  carrier  density, 
average  momentum,  and  average  velocity  in  the  ith  valley,  respectively.  The  avoaged  inter¬ 
valley  carrier  transition  rates  (v^)  are  defined  as  t 

WjVj,  =  fp^kpkj)f^k,)dk,dkj,  (5) 

where  is  the  transition  rate  from  ith  to  jth  valley.  The  intra-valley  and  inter-valley 
momentum  exchange  times  and  T„ij,  respectively)  are  defined  as 
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and  the  intra-valley  and  inter-valley  energy  exchange  times  (r^-,  r^-,  and  are  given  by 
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with 

The  energy  separation  between  jth  and  ith  valleys  is  given  by  A^=  where 

is  the  energy  gap  between  the  jth  (ith)  valley  and  valence  band,  and  hta^  is  the  inter-valley 

phonon  energy. 

Unlike  the  single-valley  model,  these  relaxation  times  or  rates  cannot  be  easily 
determined  hnom  the  phenomenological  equations  [7,8]  but  require  direct  evaluation  of  the 
integrals  in  £qs.  (S)-(7).  Instead  of  using  the  displac^-Maxwellian  distribution,  an  MKD 
[14-16]  can  be  used  to  evaluate  these  integrals  and  thus  the  relaxation  times  or  rates.  The 
concept  of  the  MKD  has  been  introduced,  in  the  previous  work  [14,16],  to  desoibe  carrier 
behavior  in  various  time  scales.  Depending  on  the  device  scale,  an  MKD  in  each  valley 
corresponding  to  the  scale  of  interest  is  chosen  to  evaluate  the  rdaxation  times.  The 
coarseness  of  the  macroscopic  description,  therefore,  depends  on  the  choice  of  the  MKD, 

The  choice  is  based  on  the  fact  that,  generally,  r>Tj>T,„  in  semiconductors,  where  r, 
r^;  and  t„  are  the  carrier  density,  energy,  and  momentum  relaxation  times,  respectively. 
For  example,  if  the  scale  of  interest  is  in  the  order  of  r^,  the  MKD  is  chosen  to  be 
= fi^t(k,n,  T)  whose  time  scale  is  described  by  t  and  r,;  thus  contains  only  information 

in  the  scale  greater  than  or  as  small  as  r^.  In  this  case,  the  relaxation  times  evaluated 
the  integrals  of  Eqs.  (5)-(7)  using  depend  only  on  the  mean  energy;  namely,  t^=tJ e  j, 
where  rqpresents  r,  r,,  or  Similarly,  in  the  scale  smaller  than  r,  (therefore  near  r,J, 
a  higher-resolution  MKD  can  be  taken  as  fi^(k)  =  and  thus 

In  this  paper,  only  the  energy-time-scale  MKD,jf^,,  is  di^ssed.  The  momentum-time-scale 
MKD,  will  be  studied  in  the  future. 

In  the  energy  time  scale,  the  time  scale  of  the  carrier  distribution  in  ith  vaUey  can  be 
changed  [15]  in  this  form: 

=  /jK*,[fcpn,(r,f),i,(r,f)]  (8) 

That  is,  in  the  energy  time  scale,  the  space-time  dependence  of  is  taken  to  result  from 
changes  in  iij,  and  Using  Eq.  (8)  in  the  BTE,  given  in  Eq.  (1),  to  change  the  scale  of 
evolution  of  the  distribution  fimction,  it  leads  to  a  ”nuum>kinetic”  equation  for  the  MKD  in 
the  scale  of  r,.  The  equation  for  the  energy-time-scale  MKD  has  been  found  to  be,  for  field 
along  z  direction 
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Eq.  (9a)  represents  a  distribution  function  that  obeys  a  steady-state  BTE  with  a  time-varying 
"equivalent*  field  of  the  ith  valley,  as  the  driving  field,  and  is  given  by  Eq.  (9b) 
as  a  function  of  Tj.  This  MKD  thus  responds  instantaneously  to  £,  -,  and  the  formation 
in  the  time  scale  smaller  than  r^,-  is  filter^.  can  be  obtained  mm  Eq.  (9a),  and  the 
relaxation  times  or  rates  in  each  vall^  can  be  evaluated  from  the  int^rals  given  in  Eqs.(S)- 
(7)  using  in  terms  of  E^^  or  c These  rates  or  times  are,  therefore,  tabulated  for 
solving  the  multi-valley  conservation  equations  given  in  Eqs.  (2)-(4). 

m.  The  ApoHcations  of  the  Developed  Transport  Model 

'  The  developed  non-equilibrium  multi-valley  transport  model  has  been  applied  to  electron 
transport  in  n-type  three- valley- GaAs  subjected  to  step  fields.  In  this  three-valley  GaAs,  a 
homogeneous  electron  density  is  assumed.  Results  obtained  from  the  multi-valley  and  sin^e- 
valley  conservation  equations  based  on  are  compared  to  those  from  the  three-valley 
Monte  Carlo  simulation. 

In  GaAs,  at  fields  below  2.0  kV/cm,  the  electron  population  in  the  upper  (L  and  20 
valleys  is  negligible,  and  thus  the  single-valley  conservation  equations  can  properiy  describe 
the  electron  transport  behavior.  The  excellent  agreement  between  average  dectron  velocities 
obtained  from  the  single-valley  conservation  equations  and  Monte  Carlo  method  is  illustrated 
in  Fig.  1  where  a  stq>  field  (from  0.2  to  2.0  kV/cm)  is  applied  to  the  n-type  GaAs  at  3(X) 
K.  However,  when  applied  ^d  increases,  the  upper-valley  population  bewmes  significant 
and  inter-valley  scattering  dominates  the  scattering  processes.  In  this  case,  the  single-valley 
description  provides  insufficient  information  on  electron  transfer  between  vaUeys  and  a  multi¬ 
valley  description  is  necessary.  The  deviation  between  average  electron  velocities  as  a 
consequrace  of  a  step  change  in  field  from  40  to  5  kV/cm  in  GaAs  determined  by  the  single¬ 
valley  description  and  Monte  Carlo  method  is  (^arly  shown  in  Fig.  2(a)  where  die  multi¬ 
valley  results  based  on  the  developed  transport  model  are  also  inclw^.  Similar  to  the 
results  in  the  previous  work  [16],  the  single-valley  results  respmid  to  the  rapid  change  in 
field  faster  than  the  Monte  C^arlo  results.  The  velocity  docribed  by  the  multi-valley 
conservation  equations  is  in  good  agreement  with  that  from  Monte  (Zado  calculatitms; 


however,  the  difference  is  still  observable.  The  difference  of  the  average  velocities  between 
the  multi-valley  conservation  equations  and  Monte  Carlo  calculations  in  Fig.  2(a)  actually 
results  from  the  deviation  of  average  velocities  of  electrons  in  F  valleys  as  illustrated  in  Fig. 
2(b).  However,  these  two  approaches  provide  nearly  the  same  evolution  of  the  L- valley 
velocity  shown  in  Fig.  2(c). 


Fig.  1  Evolution  of  average  velocity  of  electrons  in  n-typc  GaAs 
subjected  to  a  step  change  in  field. 


The  difference  of  the  F-valley  velocities  shown  in  Fig.  2(b)  determined  from  these  two 
approaches  mainly  results  from  the  oversimplified  description  for  the  macroscc^c  effective 
mass,  m*,  used  for  solving  the  multi-valley  conservation  equations.  That  is,  m  is  assumed 
to  be  a  function  of  electron  mean  energy  and  the  energy  dq)endeoce  of  m*  is  detennined  by 
^e  steady  state  relation  between  average  mommtum  and  velocity;  i.e.,  m*=^m*(7)— 
Pu(i)/u„(i),  where  e  is  the  instantaneous  mean  energy  of  electrons.  Although  this 
description  for  m*  has  been  widely  used  [8-13]  in  device  a^ysis,  it  has  been  found  in  the 
present  calculations  that  the  velocity  dependence  of  m*  might  be  significant  Fig.  2(d)  shows 
the  evolution  of  obtained  fiom  the  Monte  Cario  method  compared  to  that  of  the  energy- 
dependent  effective  mass.  It  is  noted  that  electrons  in  F  valley  described  by  mpfc)  is  lighter 
[see  Fig.  2(d)]  during  the  interval  of  velocity  undershoot  [see  Fig.  2(b)].  llie  smaller 
is,  therefore,  responsible  for  the  lar|er  velocity  [see  Fig.  2(b)]  obtained  from  the  multi-vaUey 
transport  model,  which  assumes  m  during  the  undershoot  intervaL  However,  the 

velocity  undershoot  shown  in  Fig.  2(c)  is  not  pronounced  in  L  vall^,  and  Aerefote  m^(T) 
shown  in  Fig.  2(e)  provides  an  excellent  description.  As  a  consequence,  the  L-vall^ 
velocities  detennined  by  the  multi-valley  transport  model  and  Monte  Cado  method  are  in 
excellent  agreement 


Effeclivs  Moss  (mVm.) 


Fig.  2  Evolution  of  average  velocity  and  macroscopic  effective  mass  of  electrons, 
(a)  Total  average  velocity,  (b)  T-valley  velodty,  (c)L-vallcy  velocity,  (d)  F-valley 
effective  mass,  and  (e)  L-valley  effective  mass. 
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Fig.  3  Evolution  of  average  velocity  and  macroscopic  effective  mass  of  electnms. 
(a)  Total  average  velocity,  (b)  T-valley  velocity,  (c)  L-valley  velocity,  (d)  F-valley 
effective  mass,  and  (e)  L-valley  effective  mass. 


Similar  behavior  has  also  been  observed  during  the  overshoot  interval.  The  average 
electron  velocity  in  the  n-type  GaAs  under  the  influence  of  a  step  change  in  field  from  5  to 
40  kV/cm  is  illustrated  in  Fig.  3(a)  where  three  different  approaches  are  included  for 
comparison.  Again,  the  average  velocity  described  by  the  single-valley  approximation 
reaches  steady  state  faster  than  the  multi-valley  description.  The  difference  between  the 
average  velocities  in  Fig.  3(a)  obtained  from  Ae  muld-valley  conservation  equations  and 
Monte  Carlo  calculations  results  from  the  relatively  low  overshoot  velocity  in  T  valley 
described  by  the  conservation  equations  (compared  to  the  Monte  Carlo  results),  as  shown  in 
Fig.  3(b).  The  deviation  of  F-valley  velocities  between  these  two  approaches  is,  again,  the 
consequence  of  the  oversimplified  dependence  of  ni^(7)  used  in  the  conservation  equations. 
During  the  overshoot  interval,  larger  energy-dependent  effective  mass,  [see  Fig,  3(d)] 

is  used  in  conservation  equations,  and  thus  results  in  smaller  overshoot  velocity  in  Fig.  3(b). 
However,  ml( e)  is  almost  the  same  as  the  effective  mass  obtained  from  the  Monte  Carlo 
calculations,  as  shown  in  Fig.  3(e),  since  the  velocity  overshoot  in  L  valley  is  not 
pronounced  [see  Fig.  3(c)].  This  leads  to  good  agreement  of  the  I-valley  velocities  obtained 
from  these  two  approaches. 


IV.  Conclusions 

A  non-equilibrium  multi-valley  transport  model  has  been  developed  and  used  to  study 
the  electron  transport  phenomena  in  three-valley  GaAs.  The  results  show  that  the  single¬ 
valley  transport  model  provides  an  excellent  description  for  electron  transport  in  GaAs  at 
fields  below  2.0  kV/cm  because  inter-valley  scattering  is  negligible.  At  high  fields  where 
inter-valley  ottering  is  significant,  a  model  including  information  on  carrier  transfer 
between  valleys  is  desirable.  The  developed  transport  model  can  provide  inter-valley 
transfer  effects,  as  can  be  seen  from  the  excellent  results  presented  in  Figs.  2(c)  and  3(c). 
The  results  illustrated  in  Figs.  2  and  3  have  suggested  that  the  macroscopic  efrective  mass 
might  be  strongly  dependent  on  the  average  velocity  when  velocity  overshoot  or  undershoot 
behavior  is  pronounced.  The  proper  description  for  the  macroscopic  effective  mass  will  be 
studied  in  the  future. 
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V.  Monte  Carlo  Investigation  of  Acoustic  InteivaDey  and  IntravaUey 
Deformation  Potentials  for  GaAs 


1.  Intnxiuction 

The  study  of  charge  transport  in  semiconductors  and 
devices  requires  a  number  of  material  parameters,  such  as 
the  effective  mass,  acoustic  deformation  potential,  various 
phonon  energies,  energy  separations  between  valleys,  and 
intervalley  and  intravalley  deformation  potentials,  etc. 
These  parameters  are  necessary  for  determination  of  the 
carrier  mobility,  diffusion  coefficients,  energy  and 
momentum  relaxation  times,  and  thermal  conductivity. 
However,  there  are  large  uncertainties  in  the  experimental 
and  theoretical  values  of  some  of  these  parameters.  In 
GaAs,  considerably  different  values  for  the  acoustic 
deformation  potential'  '*  in  the  F  valley  (D^),  effective 
mass'*''’  in  the  L  valley  (mj,  energy  separations'*'* 
between  F  and  L  valleys  (Atn),  and  L-L  intravalley  and  F- 
L  intervalley  deformation  potentials'*’'*  '*'’*  (Du.  and  Dn., 
respectively)  have  been  suggested  in  the  past  decade.  For 
example,  a  value  of  7.0  eV  for  of  GaAs"  '***”  has 
been  used  for  studies  of  electron  transport  in  GaAs,  or 
heterostructure  GaAs-GaAIAs.  Althwgh  this  value, 
together  with  some  adjustable  (uncertain)  material 
parameters,  has  given  reasonable  theoretical  steady-state 
results"''*^*  compared  to  the  experimental  data,”-**  recent 
experimental''*  and  theoretical”  results  have  suggested  a 
larger  value  for  D^. 

The  experiment'  of  free  electron  absorption  of  Se- 
doped  GaAs  has  given  a  value  of  1S.7  eV  for  0^.  The 
an^ysis*  of  temperature  dependence  of  mobility  in  n-type 
GaAs  has  led  to  16.0  ±0.j  eV.  Recently,  it  has  been 
found”  that  a  high-mobility  heterostnicture  exhibits  linear 
temperature  dependence  on  the  inverse  mobility  of  electrons 
at  very  low  temoerature.  This  indicates  that  the 
contribution  of  acousdc-ptionon  scattering  dominates  the 
temperature  dependence  of  mobility  at  very  low 
temperature.  The  analysis*'’  of  this  phenomena  has  yielded 
a  value  of  13.0- 14.0  eV  for  D^,  in  agreement  with  the 
13.3  eV  suggested  by  Prica.”^  Ifowever.  a  relatively 


small  value  of  D^=9.3  eV  has  been  obtained  from  the 
measurement*  of  the  band-edge  defonfiatien  potential  in 
GaAs.  More  recently,  theoretical  calculations  have  also 
arrived  at  similar  values:  8,  9.3,  and  9.6  eV  for  D^.”  The 
various  calculated  or  measured  values  for  may  be  due 
to  the  incomplete  model  and/or  experimental  errors. 
Moreover,  the  choices  of  the  values  for  the  other  material 
parameters  of  GaAs  in  F  valley  (such  as  nv  and  the  optical- 
phonon  energy  hbj  used  in  the  calculations  may  also  give 
rise  to  the  deviation  in  the  value  of  D,|.. 

However,  these  results  suggested  a  value  for 
somewhat  between  9.0  and  16.0  eV,  instead  of  7.0  eV  that 
is  usually  used.''''*^'”  Moreover,  the  change  of  the  value 
of  D,r  from  7.0  eV  will  also  necessitate  the  modification  of 
some  other  material  parameters  (including  Dn.and  D|x)  that 
were  determined"  numerically  based  on  the  Monte  Carlo 
simulations  using  D,r='^  0  eV.  In  this  paper,  Monte  Carlo 
simulations,”  using  nonparabolic  bandstnicture,  for  electron 
traiuport  in  GaAs  at  300  K  have  been  carried  out  to 
investigate  the  effect  of  the  acoustic  deformation  potential 
D,r  on  the  low-field  mobility  of  electrons,  and  thus  to  arrive 
at  a  reasonable  value  of  D^.  The  dependence  of  the 
mobility  and  velocity-field  relation  on  the  material 
parameters  including  D^-,  bb.,  mr,  m^,  Dn,,  Du,,  and  Acn, 
is  discussed  for  fields  less  than  IS.O  kV/cm.  The  influence 
of  Dj;  Dfi,  and  Du,  on  the  peak  velocity,  threshold  field 
for  Gunn  effect,  and  negative  differential  mobility  is  also 
presented.  Although  the  material  parameters  used  in  Monte 
Carlo  simulations  are  taken  to  te  constant,  it  should  be 
pointed  out  that  all  parameters  can  actually  vary  with 
energy. 

2.  Uw-Field  Mobility 

At  very  low  field,  polar  optical-phonon  scattering 
dominates  the  scattering  processes  in  GaAs  at  room 
temperature.  However,  the  low-field  mobiliQr  of  electrons, 
p,  also  strongly  depends  on  acoustic-phonon  scattering  since 
it  is  a  nearly  isotropic  scattering  in  contrast  to  the  forward 
scattering  of  polar  optical  phonons.  The  value  of  p  is, 
therefore,  strongly  dependent  on  D,r.  b  addition  to  Oj-, 
our  calculatton  alao  show  til  p  mJtibib  high  unsitivity  le 


the  value  of  optical-phonon  energy  h2>..  This  is  because  its 
value  is  comparable  to  thermal  energy  at  room  temperature, 
and  thus  an  increase  in  the  value  of  h2>.  used  in  Monte 
Carlo  calculations  at  low  fields  (where  energy  mean  energy 
is  close  to  the  thermal  energy)  will  considerably  reduce  the 
number  of  electrons  sufferi.'g  the  optical-phonon  emission. 
Moreover,  near  the  thermal  energy,  the  optical-phonon 
scattering  rate  slightly  decreases  with  the  value  of  hZ>.. 
Consequently,  the  increase  in  hS>.  can  raise  the  value  of  m- 
The  dependence  of  low-field  mobility  on  optical  phonon 
energy  is  shown  in  Fig.  1  for  different  values  of  D^.  Fig.  1 
also  illustrates  the  strong  dependence  of  low-field  mobility 
on  the  effective  mass  mr-  Results  of  ^  shown  in  Fig.l 
using  0^=9. 3  eV  and  mr=0.063m.  are  close  to  those  using 
Dj.*  13.5  eV  and  mr=0.059m..  The  value  of  D^-  can  be, 
therefore,  estimated  from  Fig.l  using  the  reliable 
experimental  results  for  effective  mass  mr,  optical-phonon 
energy  h2>.,  and  low-field  mobility 

The  experiments’^'**  of  lUI  effect  measurements 
lead  to  a  value  of  electron  Hall  mobility,  tin,  slightly  greater 
than  8600  cmVv-s  for  pure  GaAs  at  room  temperature, 
where  (based  on  the  van  der  Pauw  method*'),  and 

rn  is  the  Hall  factor  which  depends  on  magnetic  field.  The 
v^ue  of  Ch  was  measured**-**  to  be  approximately  1. 1  at  the 
magnetic  field  used  in  the  measurements.”'**  The  low-field 
mobility  of  electrons  in  GaAs  at  room  temperatxire  is, 
therefore,  calculated  to  be  approximately  7900  cmVv-s. 
This  value  is  in  good  agreement  with  the  theoretical 
values;**  **  7900  and  8000  cmVv-s.  According  to  the 
analysis**  of  the  infrared  reflection  spectra  of  GaAIAs  mixed 
crystal  using  Kiamers-Krorung  technique,  hS>.  is  found  to  be 
equal  to  36.2  meV  which  is  suggested  by  many  authors*-'*-** 
and  also  close  to  the  value  36.0  or  36.1  meV  used  in  recent 
papers.'-**-**-**-**  Besides,  the  experimental  measurements** 
have  led  to  mr/m,  >0.064  ±0.(X12  for  GaAs  at  room 
temperature,  which  is  in  good  agreement  with  the  values 
used  in  some  papen"  '*:  0.063  and  0.0638. 


Fig.  1.  Low-field  mobility  of  electrons  in  GaAs  as  a 
(unction  of  optical-phonon  energy  for  different  values  of  D,r 
and  mr.  The  effect  mass  nir  is  taken  to  be  0.063in. 
for  all  results  except  for  those  reptesenied  by  *  where 
mr»0.059m,.  Symbols:  ODj.>16.0eV;  □Dj.  =  14.0eV; 
xDrf- 13.0  eV;  0  Dj.- 12.0  eV;  ☆  Dj.- 10.0  eV; 
aDj.-9.3  eV;  ■D.r-’.OeV:  eDj.-13.3  eV. 


If  the  values  of  n,  h2>.,  and  mr  are  chosen  to  be 
8000  cm*/v-s,  36.1  meV,  and  0.063m„  respectively,  the 
acoustic  deformation  potential  D,r  of  is  to  be  as  large 
as  12.7  eV  according  to  Fig.l.  However,  the  use  of 
Dj.=  16.0  eV,  hit.=36.1  meV,  and  mr=0.063m,  gives 
M-72S0  cm*/v-s.  On  the  other  hand,  the  mobility  as  high 
as  8700  cm*/v-s  is  obtained  when  using  eV, 

h2>,=36.1  meV,  and  mr=0.063m,. 

3.  Velocity-Field  Relation  Near  and  Above  the 
Threshold  Field 

At  fields  near  the  threshold  field,  the  velocity-field 
relation  also,  in  part,  depends  on  the  value  of  D,r- 
dependence  is  illustrated  in  Fig.2.  The  material  parameters 
used  in  calculations  for  Fig.2  are  from  Refill  except  for 
D,r-  The  experimental  results**-**  are  also  included  for 
comparison.  As  shown  in  Fig.2,  in  addition  to  the  smaller 
low-field  mobility  of  electrons  n,  the  larger  D,r  >lao  leads 
to  smaller  peak  velocity,  and  larger  threshold  field.  As 
field  increases  to  a  value  above  IS.O  kV/cm,  the  calculated 
average  velocity  is  almost  independent  of  the  value  of  D,r 
because  of  much  stronger  inter^ley  scattering.  In  Fig.2, 
the  deviation  from  experimental  results  at  fields  around  the 
threshold  field  is  evid^  when  using  larger  D,r-  However, 
in  Fig.2.  a  smaller  value  of  optical-phonon  energy 
(hb,>35.36  raeV)  than  36.1  meV  (as  discutsed  in  the 
previous  section)  Ins  been  used.  Based  on  Fig.l,  the  use  of 
hb.>36.1  meV  will  raise  the  low-field  mobility,  and  thus 
will  increase  the  slope  of  the  velocity-field  curve  at  low 
field  in  Fig.2.  Also,  a  relatively  laiye  value  of  T-L 
deformation  potential  (Dn.>1.0xl(F  eV/cm)  has  been  used 
in  Fig.2.  Although  this  value  for  Dn,  hu  bm  widely  used 
in  many  papen,”-**-**'’’ values  for  Dn.  between  1.5x10*  and 
7.0x10*  eV/cm  have  been  suggested  in  recent  papen'*’*'-**'** 
(see  the  recent  review  and  theoretical  calculations  in  Refs. 
25  and  26).  In  addition,  values  for  effective  mass  m,.,  and 
the  energy  sepantion  between  F  and  L  valleys  Acn. 


Fig.  2.  Average  velocity  of  electrons  in  GaAs  as  a  finKtion 
of  electric  field  for  different  values  of  D,|.  varying  from  7.0 
eV  to  16.0  eV.  Also  included  are  the  experimental  results 
(Refe.  33  and  34). 


(0.222m,  and  0.33  eV  are  used,  respectively,  in  Fig.2)  used 
in  many  papers  are  considerably  different.  The  value  for 
mt  has  been  ranged  from  0.18m.  to  0.29m,,"  '’'””  and 
between  0.28  eV  and  0.34  eV  for 

Our  calculations  indicate  that  reducing  the  value  of 
Dn  leads  to  a  decrease  in  the  threshold  field  and  peak 
velocity.  This  is  because  electrons  are  heated  up  more 
easily  when  using  smaller  Dn..  Therefore,  energy 
distribution  of  electrons  in  the  F  valley  becomes  flatter, 
which  increases  the  number  of  electrons  with  energy  above 
A<n..  On  the  other  hand,  increasing  the  value  of  d«n. 
results  in  increases  in  the  threshold  field  and  peak  velocity. 
The  threshold  field  and  peak  velocity  are  not  sensitive  to  the 
deformation  potential,  Du,,  but  vary  slightly  with  nii.. 
Moreover,  the  negative  slope  of  the  velocity-field  curve  (ie, 
the  negative  differential  mobility)  at  fields  between  the 
threshold  field  and  approximately  8.0  kV/cm  strongly 
depends  on  DpLi  and  is  strongly  dependent  on  Du.  and  mi. 
for  fields  greater  than  7.0  kV/cm. 

Ftg.3  shows  that  the  same  threshold  field  can  be 
achieved  by  properly  varying  the  values  of  Drt  and  D,|.. 
The  use  of  larger  D^-  (13.5  eV)  is  to  result  in  a  larger 
threshold  field  which,  however,  is  compensated  by  the  use 
of  a  smaller  D^.  (4.5x10*  eV/cm).  As  a  consequence,  the 
velocity-field  curve  with  Dn.~4.5xl0'  eV/cmand  D,r~  f3.5 
eV  exhibits  the  same  threshold  as  the  one  with  Dn.~7.0xl0' 
eV/cm  and  D,r~  10.0  eV.  Based  on  the  same  reasons,  the 
use  of  Drf-~7.0  eV  and  Dn.~1.22xl0^  eV/cm  can  also 
provide  the  same  threshold.  The  peak  velocity  is  smaller 
when  using  D,r”  13.5  eV  due  to  stronger  acoustk-phonon 
scattering.  In  addition,  the  negative  slope  at  fields  4.0  and 
8.0  kV/cm  is  relatively  large  fiir  the  curve  using 
Dn.~4.5xl0*  eV/cm.  This  is  because  electrons  are  heated 
up  more  easily,  as  previously  mentioned,  and  more 
electrons  are  pumped  into  upper  (low-mobility)  valleys;  thus 
the  velocity  decreues  fester  with  applied  field. 

Fig.4  illustrates  that  Monte  Carlo  simulations  using 
A«n.~0.33  eV  and  Dn.~8xl0'  eV/cm  give  the  same  peak 
velocity  and  threshold  field  u  those  using  Aen.~0.3S  eV 


Fig.  3.  Average  veloctty  of  electrons  in  QaAs  as  a  ftmction 
of  electric  field  for  different  values  of  D,r  and  Dn..  Also 
included  are  die  experimental  results  (lUfe.  33  and  34). 


Fig.  4.  Average  velocity  of  electrons  in  GaAs  as  a  function 
of  electric  field  for  diffinem  values  of  Acn,  Dn,  and  Du.. 
Also  included  are  the  experimental  results  (Refe.  33  and 
34). 


and  Dn~d>0xl0'  eV/cm.  Same  values  for  mr,  mt,  D,r> 
and  hb,  are  used  in  the  simulations  for  Fig.4.  The  use  of 
larger  Dn  (8.0x10*  eV/cm)  intends  to  increase  both  of  the 
peak  velocity  and  the  threshold  which  can  be  balanced  by 
using  a  ampler  value  of  Asn  (0>33  eV).  Again,  the 
negative  slope  at  fields  between  4.0  and  7.5  kV/cm  is  larger 
for  the  smaller  Dn>  The  dependence  of  the  negative  slope 
on  the  value  of  Du.  at  fields  greater  than  7.0  kV/cm  is  also 
illustrated  in  Fig.4.  The  average  velocity  with  the  value  of 
Du.-8  .0x10*  eV/cm  decreases  fester  than  the  one  with 
Du.~6x10*  eV/cm  as  field  increases;  both  uses  the  same 
Acn  and  Dn-  This  is  because  the  larger  Du.  causes  higher 
scattering  rate  in  the  L  valleys. 

4.  Conclusions 

In  this  paper,  the  acoustic  defermation  potential, 
D,P,  has  ben  estimated  to  be  12.7  eV  according  to  the 
dependence  of  the  low-field  mobility  on  D,r  and  hb,,  as 
shown  in  Fig.l.  However,  the  information  provided  by  the 
available  experimental  measurements  is  not  sufficient  to 
determine  intervalley  and  intravalley  deformation 
potentials.  As  can  be  seen  in  Fig.3  and  4,  proper 
adjustments  of  the  uncertain  parameters  (nii.,  Actl,  Do,  and 
Du)  can  provide  reasonable  Monte  (ferlo  results  of  the 
velocity-field  curve  to  fit  the  experimental  results.”-” 
There  is  no  unique  set  of  the  parameSers  that  can  generate 
results  in  good  agreement  with  experimental  measurements 
of  the  velocity-field  curve.  In  order  to  ftirther  restrict  the 
set,  it  it,  therefore,  desirable  to  also  include  the  comparison 
between  experimental  and  theoretical  results  of  the 
longiudinal  diffioioo  coefficient  which  strongly  depends  on 

the  value  of  mi.  Acru  Da  *0*1  ^ 

the  threshold  field.  However,  the  available  experimental 
results**-**-”  for  the  tfiffiision  coefficient  ate  evidently 
different  For  example,  the  longitudinal  diffusion 
coefficieat  was  suggested  to  be  above  800.0  cm*/sec  at  3.0 
kV/cm  in  Ref.33  and  less  than  450  ctn*/tec  in  Refe.  53  and 


54.  Moreover,  as  previously  mentioned,  the  experimental 
and  theoretical  values  for  the  effective  mass  mi^,  and  the 
energy  separation  Atn.  that  have  been  used  in  previous 
studies,  are  considerably  different.  To  determine  the  value 
for  Dn.  and  Du.,  more  reliable  experimental  results  for 
longitudinal  diffusion  coefficient,  the  energy  separation 
between  F  and  L  valleys,  and  the  L-valley  effective  mass 
are  desirable. 

In  submicron  and  high-speed  semiconductor  devices, 
carrier  transport  behavior  is  characterized  by  energy  and 
momentum  relaxation  times’^”  which  strongly  depend  on 


intervalley  and  intravalley  scattering  at  fields  above  the 
threshold,  and  thus  on  the  values  of  Dn.  and  Du..  To 
develop  a  transport  model  for  these  devices,  and  to  study 
carrier  dynamics  in  &st  transient  situations,  it  is  necessary 
to  obtain  more  accurate  values  for  Dn.  and  Du.. 
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1.0  INTRODUCTION 

Electro-optic  sampling  [Valdmanis  et  al  1983,  Kolner  and  Bloom,  1986] 
measures  electric  field  magnitudes  with  both  high  spatial  and  temporal 
resolution  via  a  medium  with  electric  field  dependent  refractive  indices. 
These  characteristics  provide  a  first  opportunity  to  experimentally 
investigate  short  time-scale  phenomena,  including  velocity  overshoot  [Nuss 
et  al.,  1987,  Meyer  et  al.  1988],  local  field  vamiations  in  semiconductor 
materials,  and  opto-electronic  switching  not  previously  observable  directly. 
The  technique  has  been  applied  to  the  observation  of  pulse  propagation  in 
GaiAs  integrated  circuits  where  the  GclAs  is  used  as  both  the  circuit 

substrate  and  the  electro-optic  medium  [Weisenfeld  et  al.  1987,  Zhang  and 
Jain  (1987),  Rodwell  et  al.  1986],  and  subpicosecond  photoconductivity 

(Iverson  et  al.,  1988,  Chamoun  et  al.,  19891.  Space-charge  phenomena  are  key 
problems  to  the  understanding  of  interf  acial  polarization  and  electret 

formation,  carrier  transport  across  metal-semiconductor,  and  metal-insulator 
barriers,  transition  from  electrode  limited  to  bulk  limited  conduction 
processes  in  metal-semiconductor-metal  systems,  and  to  the  imderstanding  of 
the  self-sustained  conduction,  known  as  lock-on  [Brinkmann  et  ad  1991],  seen 
in  semiconductors  at  high  fields.  Non-equilibrium  charge  carrier  tramsport 
in  semiconductors  yields  information  about  camrier-phonon  (Kutt  et  al,  1991] 
and  camrier-cairrier  interactions  and  high  field  phenomena  like  velocity 

overshoot,  ballistic  transport,  and  resonamt  tunneling.  Previously  used 
techniques  have  been  either  electronic  or  optical  but  but  lack  the  time 
resolutions  adequate  to  observe  relaxation  phenomena  on  the  tens  of 
femtoseconds  to  picosecond  time-scales. 

2.0  THEORETICAL  AND  DESIGN  ASPECTS 

An  apparatus  for  measurement  of  the  field  strengths  of  short  pulses 
produced  by  photocarriers  in  semiconductors  requires,  for  many  applications, 
sensitivity  to  small  fields  and  high  spatial  and  temporal  resolution.  These 
characteristics  can  be  obtauned  from  a  combination  of  the  high  frequency 
response  of  the  polarizability  of  electro-optic  materials  (such  as  LiTaO^ 
and  LiNbO^),  ultra-short  pulses  of  modelocked  laser  systems  and  narrow 
bandpass  detection  systems.  This  section  describes  the  principles  upon  which 
electro-optic  sampling  is  based  and  methods  available  to  implement  them. 
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2.1  The  Index  Ellipsoid  and  Impermeability  Tensor 

The  electric  field  at  a  point  in  a  dielectric  is  the  sum  of  two  fields: 
D/e  associated  with  free  charges  and  -P/c  associated  with  bound  charges, 

O  O 

that  is, 

D  =  e  E  +  P.  (1) 

O 

where  P,  the  electric  polarization  is  the  dipole  moment  per  unit  volume  at  a 
given  point  and  is  a  macroscopic  quantity.  In  general,  P  is  a  tensor  and 
can  be  expressed  as 

P  =  e(yE  +  >E  ■•■>;£), 

X  o  '^11  X  '*■12  y  '*■13  z  ’ 

P  =  c  (z  E  +  y  E  +  X  E  (21 

y  o  '^21  X  ^22  y  *23  z 

P  -  e  (;t  E  +  y  ,E  +  y  ,E  ) 

z  o  31  X  32  y  33  z 


allowing  D  to  be  expressed  in  terms  of  £: 
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or  D  »  c  E  where  c  *  c  (1  +  x  )  and  for  c  a  real  dielectric  tensor 

1  ij  j  ij  o  ij 

c  =  e  . 

U  Ji 

The  energy  density  in  the  stored  electric  field  in  an  anisotropic 
medium  cam  be  written  in  terms  of  the  electric  field  vector  and  the 

displacement  tensor 

U  =  i  E.D.  (4) 

«  2 

From  this  the  surfaces  of  constant  energy  in  D  space  are  determined: 


«  2U  , 

0 


where  c  ,c  ,  and  e  are  the  principal  dielectric  constants.  If  we  replace 

X  y  z 

D/2U  by  r  and  define  the  principle  indices  of  refraction  n  ,  n  ,  and  n  by 
«  xyz 

2 

previous  equation  becomes 


-2  2  2 
n  n  n 

xyz 

This  is  the  equation  of  an  ellipsoid,  known  as  the  index  ellipsoid  or 
optical  indicatrix.  It  determines  the  two  indices  of  refraction  amd  the  two 
corresponding  directions  of  D  associated  with  the  two  independent  plane 
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waves  that  ceui  propagate  along  an  arbitrary  direction  s  in  the  crystal. 

The  impermeability  tensor  is  related  to  the  refractive  index  tensor  by 
7)^^  =  e^(e  and  the  quantum  theory  of  solids  predicts  that  the  optical 
dielectric  impermeability  tensor  depends  on  the  distribution  of  charges  in 
the  crystal.  Hence,  the  application  of  an  electric  field  will  result 
in  a 

redistribution  of  bond  charges  and  possibly  a  slight  deformation  of  the  ion 
lattice,  resulting  in  a  change  in  the  optical  impermeability  tensor.  The 
phenomena  is  known  as  the  electro-optic  effect  and  the  electro-optic 
coefficients  are  traditionally  defined  by 

where  E  is  the  applied  electric  field,  the  constants  r^^  are  the  linear  (or 
Pockels)  electro-optic  coefficients,  and  the  s^^^  are  the  quadratic  (or 

Kerr)  electro-optic  coefficients. 

In  the  principle  axis  system  for  the  impermeability  tensor,  c  is 

a  diagonal  matrix.  So 
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(8) 


On  applying  an  electric  field  to  a  non-centrosymmetric  crjrstal  the  quadratic 
term  is  assumed  to  of  much  smaller  magnitude  than  the  linear  term  and  the 
index  ellipsoid  becomes 
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where  5  is  the  Kronecker  delta: 
ij 
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For  a  lossless  and  optically  inactive  medium  the  dielectric  tensor 
is  symmetric  and  so  the  impermeability  tensor  must  also  be  symmetric. 
Consequently  the  indices  i  zmd  j  can  be  permuted  in  the  above  equation. 
The  permutation  symmetries  are  summarized  ais 

r  =  r 

IJk  Jlk 

Because  of  the  permutation  sjrmmetries,  it  is  convenient  to  introduce 
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contracted  indices  to 
where 


abbreviate  the  notation.  Using  this  notation  r  =>  r 

ijk  mk 


ij  11  22  33  23  13  12 

ml  2  3  4  5  6. 


LiTaO^  is  of  the  rhombohedral  crystal  class,  R3C  or  space  group,  3m 
or  point  group  and  the  electro-optic  tensor  has  the  form  [Yariv  and  Yeh, 
19841 


r 

mk 


LlTaO 
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(10) 


The  contraction  of  the  indices  is  a  matter  of  convenience  and  the 
matrix  elements  (6x3)  do  not  have  the  usual  tensor  transformation  or 
multiplication  properties. 

When  an  electric  field  is  applied,  the  crystal  equation  (9)  will,  in 
general,  no  longer  represent  the  principle  axis  system.  In  a  t3rpical 
alignment,  where  the  field  is  E^,  the  index  ellipsoid  is  described  by 
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For  LiTaO^  this  expression  reduces  to 
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The  effective  refractive  indices  are  field  dependent,  providing  a 
means  of  altering  the  relative  phase  of  the  E  field  components  of 
beam  propagating  perpendicular  to  the  z  direction. 


(12) 

practical 
the  light 
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2.2  Transverse  Field  Light  Modulators 

The  phase  change  between  polarizations  can  be  determined  by  placing 
the  crystal  between  two  linear  polarizers  in  a  configuration  known  as  a 
Pockets  cell  as  shown  in  Fig.l  .  A  linearly  polarized  beam  of  intensity 


Fig.l  A  Pockel’s  cell.  The  relative  phase  shift  introduced  to  the  two 
perpendicular  polarizations  by  the  field  dependent  birefringence  alters  the 
intensity  I  transmitted  through  polarizer  P^. 

enters  the  crystal  and  after  propagating  a  length  L  through  the  applied 
field  has  a  relative  phase  shift  P  between  the  two  polarizations  aligned 
with  the  principle  axes.  The  output  intensity  of  the  analyzer  can  be 
determined  using  Jones  calculus,  in  which  matrices  represent  the  amplitude 
and  phase  changes  introduced  to  the  beam  by  individual  optical  components. 
For  the  case  of  a  phase  inducing  material  between  two  crossed  polarizers 
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Hence  the  intensity  of  the  transmitted  beam  is 

I  «  I  sin^(r/2) 

0 

where  P  *  ^  ^  »  *^nL(n’  -  n')  and  k  is  a  constant.  For  LiTaO, 


1/n’^  ■  1/n*  ♦  r  E  ,  and  l/n’*  «  1/n^ 

X  X  13  Z  Z  X 
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From  this,  an  expression  for  the  field  dependent  refractive  index  can  be 


obtained; 
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So  the  phase  shift  introduced  between  the  two  polarizations  of  a  beam 
travelling  in  the  y  direction  is 

r  =  |lE(n  -  n  )L  -  ^n^r  -  nV  )E  L  (15) 

A  e  o  Ae33  ol3z 

where  L  is  the  distance  travelled  through  the  field  and  A  is  the  free  space 
wavelength. 

2.3  Half-wave  Voltage 

The  half  wave  volteige,  V  is  a  useful  parameter  for  characterizing 
electro-optic  modulators.  It  is  the  voltage  required  to  introduce  a  phase 

shift  between  the  emerging  orthogonal  polarizations  by  n  radians  and  is 

often  quoted  as  a  measure  of  Pocicels  cell  responsivity. 

V  *  Ad/(L(n^r  -n^r  ))  (16) 

a  •  33  o  13 

where  d  is  the  separation  of  the  electrodes.  For  the  LiTaO^  modulator 

described  in  section  3,  with  d  =  0.7  mm,  L  taken  to  be  0.3  mm,  A  *  580  nm 
and  n^r  -  n^r  =  2.285  e-10  m/V  the  half-wave  voltage  is 

9  33  o  13 

V  =  5.9  KV. 
n 

2.4  Quarter-wave  Phase-shift  Bias 

From  equ  14  the  intensity  change  as  a  function  of  phase  shift  is  maximum 
when  the  F  is  a/2  and  at  this  phase  shift  the  relationship  is  approximately 
linear  as  shown  in  Fig. 2.  Advantage  of  these  characteristics  can  be  made  by 
introducing  an  adjustable  (Soliel-Babinet)  compiensator  into  the  beam  path 
and  tuning  to  introduce  a  constant  phase  shift  »  a/2.  In  this  way  the 
sensitivity  to  small  phase  shifts,  AF  =  aV/V^,  introduced  by  variations  in 
field  strength,  given  by 

I  =  I  sin^((F  +AF)/2)  (17) 

0  0 

is  maximized  and  most  linear.  Here  V  is  the  voltage  associated  with  the 
field  E  .  Note  that  F  *  a/2  and  small  AF,  I  “  (I  /2)(1+AF). 

z  0  0 
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1.20  -a 


Fig.  2  A  linear  response  the  transmitted  beam  is  obtained  for  small  phase 
changes  at  the  quarter  wave  bias  point  (F^  =  it/2). 

2.S  Mlcrostrip  Design 

The  field  E  need  not  be  constant.  In  fact  the  rapid  time  variaticm  of 
z 

the  field  is  the  motivation  for  the  use  of  electro-optic  sampling.  The 

response  time  of  the  polarization  of  the  crystad  is  of  the  order  of 

femtoseconds,  hence  the  modulation  can  occur  on  short  time  scales,  providing 

an  opportunity  for  characterization  of  rapidly  changing  fields. 

To  provide  an  accurate  means  of  determining  short  duration  field 

pulses,  propagation  from  the  point  of  origin  to  the  probing  location  within 

the  LiTaO  must  occur  with  minimum  distortion  of  the  field.  Such  a 
3 

requirement  places  constraints  on  the  wave-guiding  structure  used.  In 
general,  wave  guides  support  more  than  one  mode,  and  the  phase  velocity  of 
each  mode  is  frequency  dependent  giving  rise  to  dispersion.  Any 
discontinuities  in  the  wave-guide  excite  all  modes  and  give  rise  to 
reflections.  Of  the  numerous  waveguiding  structures  in  use,  balanced 
stripline  (Fig.3)  and  coplanar  stripline  provide  suitable  transmission  line 
geometries.  The  case  of  a  balanced  stripline  is  described  here.  Because  of 
its  symmetry  it  can  be  analyzed  as  a  mlcrostrip  of  height  h  in  which  the 
LiTaO  serves  as  a  substrate  for  thin  film  metal  conductors  of  width  w. 

3 

Since  the  field  lines  between  the  narrow  conducting  strip  and  the  ground 
plane  do  not  lie  entirely  inside  the  substrate  the  propagating  mode  is  not 
truly  transverse  electromagnetic  (TEIM)  but  quasi-TEM.  The  field  lines  lie  in 
two  media,  the  substrate  and  the  vacuum,  and  as  a  result  the  relative 
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dielectric  constant,  e  ,  does  not  provide  an  accurate  description  of  the 

r 

phase  velocity  and  is  replaced  by  the  effective  dielectric  constamt,  c 
which  is  smaller  than  e  . 


“I 

Fig.3  (a)  Balanced  stripline  (b)  llicroetrip. 

The  phase  velocity  in  the  microstrip  is  given  by 

V  =»  c/c 

p  aff 

and  the  transmission  line  characteristic  impedance  is 

Z  =  1/v  C 

0  p 

where  C  is  the  capacitance  per  unit  length  of  the  line. 

In  the  first  electro-optic  measurements  made  the  pulse  source  was 
external  to  the  microstrip  and  standard  50  Q  connectors  and  coaxial  lines 
were  used  to  connect  the  two.  To  minimize  the  discontinuity  at  the 
transition.  The  distance  between  the  inner  and  outer  conductors  of  the 
connector  was  be  matched  to  the  height  of  the  stripline,  2h. 

Closed  forms  of  expressions  for  the  determination  of  transmission  line 
parameters  have  been  determined  by  a  number  of  authors  and  are  outlined  in 
the  review  of  Bahl  and  Trivedi  (1977).  Hammerstad  (1975)  provides  an 
accurate  expression  for  determination  of  the  characteristic  impedance  in 
terms  of  the  relative  dielectric  constant. 

For  w/h  s  1, 

»  (eO/yr^^)  ln(8.h/w  +  0.25  W/h) 

where 

e  +1  c  -1 

e  ,  =  ((l+12.h/w)'‘^+0.04(l-w/h)^l 

eff  2  2 

and  for  w/h  s  1, 
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where 


Z 

0 


120  ff/v  eff 


w/h  +1.393  +  0.667  ln(w/h  +  1.444) 


e  +1  e  -1 

c  *  s —  (l+12.h/w) 

eff  2  2 


(18) 


The  w/h  ratio  obtained  from  Hammerstad’s  expression  are  plotted  in 

Fig.  4  as  a  function  of  Z  . 

0 


Fig.  4  (a)  The  calculated  microstrip  impedance,  amd  (b)  the  effective 
dielectric  constant  as  a  function  of  the  microstrip  height  to  width  ratio  as 
predicted  by  Hammerstad.  1975. 


At  higher  frequencies  the  electric  field  lines  are  contained  more 

within  the  substrate  and  e  tends  to  c  which  is  frequency  dependent.  For 

•ff  r 

waveforms  with  a  broad  frequency  spectrum  the  varying  phase  velocities 
result  in  a  change  in  the  pulse  shape  with  the  distance  travelled. 
Maintaining  short  distances  between  the  pulse  source  and  the  probe  location 
are,  therefore.  Important  to  minimize  the  effects  of  dispersion. 


2.6  Velocity  Matching  of  the  Pulse  and  Probe  Beam. 

As  described  so  far  the  pulse  and  probe  beeun  propagate  perpendicular  to 
one  another  in  the  crystal.  In  this  geometry  the  length  of  the  field  pulse 
interacting  with  the  probe  beam  is  not  simply  the  beam  diameter,  account 
must  also  be  taken  of  the  propagation  of  the  pulse  as  it  is  traversed  by  the 
probe.  This  is  illustrated  in  Fig. 5  where  the  probe  enters  the  pulse  at  A 
when  the  field  location  is  specified  by  E^.  At  time  t  later,  t  =  d/o’ 
=  (W/cosa)(n/c)  where  d  is  the  length  travelled  by  the  field  pulse  at 
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velocity  «'  in  the  crystal  of  refractive  index  n  by  the  probe,  the  probe 
emerges  at  position  B  and  the  field  has  propagated  to  E^.  The  field  pulse 
travels  a  distance  d^  at  speed  m  in  time  t,  and  the  relative  pulse  and  probe 
speed  in  the  x  direction  is  Therefore  a  probe  pulse  of  a  zero  width 


beam  integrates  features  extending  over  a  time  x. 


1/2 

in  the  field.  Substituting  «’  =  (n/clsina  and  «  =  we 

obtain 


T  =  (n  -  v'e  '  sina)  -  (20) 

•  ff  c.  cosa 

where  c  is  the  vacuum  velocity  of  light.  For  the  case  of  zero  integration, 
i.e.  T  =  0 


n  =  '  sina.  (21) 

6f  f 

Snell’s  Law  determines  the  angle  of  incidence  of  the  probe  beam: 

sin^  =  nV/e^^’^.  (22) 

For  the  average  value  of  n  and  n  =  2.1856  and  c  »  26.0,  cadculated  from 

o  «  eff 

equ  18,  8  =  69.5  degrees.  The  case  where  a  *  0  leads  to  x  *  nW/c  «  2.1856 
ps,  indicating  that  for  pulse  widths  <  20  p«  velocity  matching  will 
significantly  improve  the  time  resolution. 


Fig.  5  Velocity  matching  of  the  field  and  probe  pulses  by  adjustment  of  the 
angle  8  improves  temporal  resolution. 


2.7  Characterization  and  Reduction  of  Noise. 

Noise  on  a  signal  is  usually  detrimental,  imposing  a  practical  limit 
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to  the  performance  of  the  circuit.  Noise  is  due  to  inherent  effects  in  the 
system  rather  than  spurious  effects  such  as  a  faulty  contact  which  could  be 
eliminated.  The  three  most  commonly  encountered  types  of  noise  are  thermal, 
shot  amd  1/f  noise  [Buckingham,  1983). 

Thermal  noise  arises  f  rom  *he  random  velocity  fluctuations  of  the 
charge  carriers  (electrons  and  holes)  in  the  resistive  material.  The 
mecheuiism  is  Brownian  motion  of  the  charge  carriers  due  to  the  thermal 
energy  of  the  material.  Thermal  noise  is  present  when  the  resistive  element 
is  in  thermal  equilibrium  with  its  surroundings  and  is  often  referred  to  as 
Johnson  noise.  The  open  circuit  noise  voltage  of  a  resistor  has  a  flat 
spectrum,  that  is  the  noise  power  in  each  hertz  of  frequency  is  a  constant. 
The  noise  generated  by  a  resistance  R^^at  temperature  T  is 

V_^^(rms)  =  (4icTR|^B)’^^  (23) 

or  written  in  terms  of  the  mean  square  current 

=  4JCTB/R,  ,  (24) 

JN  L 

where  k  is  Boltzmann’s  constant  and  B  is  the  bandwidth  of  the  detector  in 
hertz.  (At  20°C  4kT  =  1.62  x  10'“v^/Hz.Q) 

Shot  noise  is  associated  with  the  passage  of  current  across  a  barrier, 
and  in  this  sense  is  a  non-equilibrium  form  of  noise.  Shot  noise  is 

frequently  encovmtered  in  solid-state  devices,  whenever  a  net  current  flows 
across  a  potential  barrier  such  as  a  depletion  layer  of  a  p-n  junction,  to 
illustrate  the  nature  of  shot  noise,  consider  the  example  of  a  thermionic 
diode  in  which  the  electrons  are  randomly  emitted  from  the  cathode  and  then 
drift  across  to  the  anode  under  the  influence  of  an  electric  field.  The 

current  associated  with  the  stream  of  electrons  fluctuates  randomly  about  a 
mean  value.  If  the  charges  act  independently  of  one  another,  then  the 

fluctuating  current  is  given  by 

‘sH  =  (25) 

where  q  is  the  electron  charge. 

1/f  noise  is  a  noise  found  in  all  systems  and  has  been  resilient  to 

theoretical  understamding  but  must  be  considered  when  establishing  the  noise 
level  of  a  system.  It  is  so  named  because  the  power  density  spectrum  vamies 
with  frequency  f  as  |f|  ,  where  a  lies  in  the  range  0.8  <  a  <  1.2. 

All  forms  of  noise  are  present  to  some  degree  at  all  frequencies. 
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hov  ^ver,  at  high  enough  detection  frequencies  the  contribution  by  1/f  noise 
is  negligible.  The  other  contributions  to  noise  can  be  reduced  by  narrow 
bandwidth  measuring  techniques.  Assuming,  temporarily,  that  the  1/f  noise 
can  be  neglected  the  minimum  detectable  voltage  can  be  determined  from  the 
signal  to  noise  ratic  where  noise  contributions  are  of  the  Johnson  and  shot 
types.  For  a  photodiode  current  in  a  load  resistor  the  signal  power  is 


(26) 

(27) 


S/N  =  R  (28) 

_ «lt  L _ 

4)cTB  +  2qt  BR 

^  DC  L 

An  illustrative  measure  of  the  signal  to  noise  ratio  (S/N)  is  obtained  by 
considering  a  particular  form  of  signal,  for  example,  the  case  of  the 
transmission  line  driven  by  a  signal  generator  of  output  V^sin  ci>t  of  small 
amplitude. 

The  modulated  probe  beam  intensity  of  equ  17  can  be  written  as 

I  =  (I /2)(l-cos(r  +Ar))  (29) 

0  0 

which  for  small  Af  reduce-  *o 

I  =(I /2)(l-cosr  +Arsinr  ).  (30) 

0  0  0 

When  the  transmission  line  is  driven  by  a  signal  generator  of  output 

V  sin  (i>  t 
0  m 

AF  *  (wV  /V  )  sin  u  t,  (31) 

ox  m 

hence  the  total  current  produced  by  a  detection  photodiode  (see  section  3) 

is 

I  xV 

I  (1-cos  r  +  sin  «  t  sin  F^l.  (32) 

tot  2  0  V  m  0 

X 

The  average  and  time  varying  signals  can  be  separated  to  yield 

I 

rj  -  I  sin^(r  /2)  (33) 

DC  2  0  0 
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sin  u  t  sin  f  . 

m  0 


(34) 


i  itV 
_  0  0 

'■.1*  "2  V 


The  mean  squared  shot-noise  current  density  is  thus 


=  2qi  =  2qt  [singer /2)] 

SN  DC  0  0 


and  the  mean  squeu-ed  signal  current  is 


„  2 

—  £  irV 

j2  0  0  .  2,_ 

‘.1.  “  r  r- 

IT 


(35) 


(36) 


«2  ,,  2 
i  irV 

8^^  sin^W2) 

^  S/N  =  - -  (37) 

(2q£  sin  (r  /2)  +  4XT/R  )B 

0  0  L 

Setting  S/N  =  1, 


0  min 


4V 

n 

f(t  sin  (F  /2)  +  2JtT/qR  )qBl 

0  0  L 

ni 

0 

sinV 

^  0  > 

(38) 


For  r^=  ii/2,  1  Kfl,  at  300  K  2kT/qRj^  =  .51e-4  and  the  expression  is 

dominated  by  the  first  term  in  the  numerator,  so 


(V  )  = 

0  min  n 


n  2 


i/qB/t 

^  DC 


volts. 


(39) 


For  a  detection  bandwidth  of  B  a  l  Hz,  determined  from  the  responsivity 
of  the  photodiode,  and  the  of  section  2.3  the  minimum  detectable  signal 
in  the  absence  of  1/f  noise  is  (V  )  *  80  uV. 

0  min 


2.8  Signal  Measurement. 

Probe  beam  intensities,  I,  are  determined  from  the  photocurrent  response 
of  a  photodiode.  The  voltage  developed  across  the  load  resistor  R^  in 
response  to  the  photocurrent  is  measured  using  the  phase  sensitive  detection 
of  a  lockin  amplifier.  This  method  increases  the  signal  to  noise  ratio  via 
bandwidth  narrowing.  The  signal  to  be  measured  is  superimposed  on  noise  and 
a  means  of  preserving  the  signal  while  reducing  the  amount  of  (broadband) 
noise  accepted  is  required.  The  use  of  a  low-pass  filter  may  work  in  some 
cases  but  where  signal  to  noise  ratios  are  low  this  is  inadequate.  A 
combination  of  modulation,  to  shift  the  signal  to  a  low  density  part  of  the 
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noise  spectrum,  and  signal  averaging,  to  reduce  the  bandwidth,  provides  an 
increased  signal  to  noise  ratio.  Bandwidth  reduction  is  accomplished  by  long 
periods  of  data  accumulation,  reducing  the  effective  bandwidth  to  Af=l/T, 
where  t  is  the  time  constant. 

Filtering  the  shifted  signal  with  a  bandpass  filter  of  high  Q  would 
require  accurate  tuning  of  the  center  frequency  and  a  slight  mistuning  would 
result  in  large  errors.  To  overcome  this  problem  the  lockin  amplifier  first 
rectifies  the  incoming  signal  and  then  filters  it.  To  accurately  rectify,  a 
reference  signal  synchronized  to  the  modulation  signal  is  required.  The 
phase  detector  circuit  of  the  lockin  amplifier  is  shown  schematically  in 
Fig.  6. 

10  k 


Fig.  6  The  reference  signal  to  the  FET  switch  rectifies  the  lockin  amplifier 
input  signal  prior  to  low  pass  filtering  that  enhances  the  signal  to  noise 
ratio. 


Assume  that  the  applied  signal  is  E^cosiut-i-^)  and  that  the  reference 

signal  is  a  square  wave  with  transitions  at  the  zeros  of  sin(wt). 

Recitification  is  achieved  by  switching  and  the  output  is  averaged  by  a 
low-paiss  filter  of  time  constant  greater  than  one  period: 

T  =  RC  »  T  =  2ii/« 

Then  the  low-pass  filter  output  is 

-2E 

<E  cos(ut+^)>|"  -  <E^cos(ut+^)>|^*J^  =  - -  sin^.  (40) 

w 

The  phase  difference  0  is  between  the  input  and  the  reference  signals.  The 
argument  extends  to  more  general  periodic  input  waveforms.  Hence  the 

bandwidth  narrowing  capabilities  of  the  lockin  amplifier  significantly 

increase  the  signal  to  noise  ratio. 
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3.0  EXPERIMENTAL  APPARATUS  AND  MEASUREMENTS. 


Fig.  7  illustrates  the  arrangement  of  components  used  for  electro-optic 
sampling.  The  synchronous  pump  and  probe  beams  originate  from  a  cw 
modelocked  (82  MHz)  Ar  ion  (Spectra-Physics  model  171-07)  pumped  dye  laser 
(Spectra-Physics  model  375).  The  dye  type  was  rhodamine  6G  and  output  pulses 
were  of  a  7ps  FWHM  autocorrelated  and  the  wavelength  range  was  560-600  nm. 
Glam-laser  double  escape  window  polarizing  prisms  (Karl  Launbrecht)  selected 
for  extinction  ratio  of  le-6  were  oriented  at  90  degrees  to  one  another  and 
functioned  as  polarizer  and  analyzer  of  the  Pockel’s  cell.  These 

prisms  transmit  one  polarization  and  reflect  the  second  at  approximately  65 
degrees  to  the  opticad  axis,  rendering  both  polarizations  avaulable  for 
determination  of  the  phase  change  of  the  beam  (see  appendix  A). 

Microscope  eyepiece  lenses  of  focal  length  11  mm  focused  the  beam  into 
the  LiTaO^.  The  free  space  beam  waist  can  be  determined  from  the  results  of 
paraxial  theory  applied  to  Gaussian  beams  (Kogelnick,  1965)  and  is  well 
approximated  by 


w 


w 


f/z 


“(I  ♦ 


(41) 


2 

where  z  irw  /A,  w  is  the  beam  diameter  incident  on  the  focusing  lens  and  X 
0  0 

is  the  wavelength.  The  beam  diameter  was  determined  using  the  knife-edge 
techniques  of  Khosrofian  and  Garetz,  1983  and  Suzaki  and  Tachibana,  1975 
which  infer  the  diameter  from  the  intensity  propagating  beyond  a  partially 
blocking  knife  edge  as  a  function  of  knife-edge  position.  Respectively  the 
beam  diameters  obtained  were  0.98  and  1.01  mm.  Substituting  the  average  of 
these  into  equ  41  the  beam  waist  diameter  obtained,  Zw,  is  4.4  fim. 

Aluminum  electrodes  of  width  0.3  mm  were  evaporatively  deposited  on  the 
LiTaO^  to  form  the  balanced  transmission  line  with  the  large  component  of 
the  field  in  the  direction  of  the  optical  axis  of  the  crystal.  The  substrate 
dimensions  were  0.5  x  0.7  X  10.0  mm  and  the  two  narrow  faces  were  optically 
polished  for  transparency  to  the  probe  beam.  The  remaining  faces  were  lapped 
flat.  LiTaO  .was  chosen  in  preference  to  LiNbO  because  of  its  larger 
electro-optic  response  to  a  given  applied  field.  Fig.  8(a).  llie  substrate 
was  mounted  with  one  conducting  film  in  contact  with  a  brass  support 
electrode  of  the  same  width.  At  each  end  a  27  GHz  SMA  connector  was  used  to 
connect  to  the  applied  bias  and  a  terminating  resistor  as  shown  in  Fig.  9. 
Termination  and  signal  display  were  performed  by  a  high  speed  oscilloscope. 


16 


£  e 

3  c  o  w  ®> 
•*«•-«« 

a>E*<=_ 

«  o  e  a  ^ 

®  u  £  E  £ 

®  •.  a  5  tt  ^ 

0  £  o  e  o>  a 

—  ■a  T  ••  £  o 

3  e  A  e  s  £  , 

«  ?  z  £  5  “ 

si's  s  I 
“•■5  3=  £  S 
^  5  i  « 

^  s  ‘o  ^  "w  s ' 


ELECTRO-OPTIC  SAMPLING  SYSTEM 


(a)  (b) 


Fig. 8  (a)  Field  dependent  birefringence  and  (b)static  bir^ringence  (b)  of 

LiTaO  2uid  LiNbO  . 

3  3 


Fig.  9  Connector  and  mount  geometry  of  the  eleetro-epUc  (LiTaO^) 
transmission  line. 

Both  LiNbO  amd  LiTaO  have  non-zero  birefringence  when  there  is  zero 

3  3 

applied  field,  i.e.  they  have  non-zero  static  birefringence,  fig  8{b).  A 
Soliel-Babinet  compensator  was  used  to  both  zero  the  phase  shift  introduced 
by  the  birefringence  and  to  set  the  phase  shift  to  the  quainter  wave  bias 
point  described  in  section  2.4. 

FND-100  (EC&G  Photon  Devices)  photodiodes  were  selected  to  detect  the 
emerging  polarization  intensities  because  of  there  large  surface  area  and 
high  responsivity  in  the  visible  part  of  the  spectrum.  They  were  reverse 
biased  (90  volts),  and  the  voltage  generated  across  a  1  Kfl  load  resistor  was 
monitored  on  a  DC  voltmeter  (Hewlett  Packard  EI2373A)  to  determine  the  phase 
bias,  see  Fig.  10.  The  output  was  also  AC  coupled  with  a  0.1  pF  capacitor  and 
the  periodic  waveform  amplitude,  due  to  the  polarization  change,  measured 
using  a  lockin  amplifier  (Standford  Research,  SR510) 
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C-0.1  uF 


FND-100 
90  V 


To  Lockin  Amplifier 


V  -  DC  volimeter 

Fig.  10  Photodiode  bias  circuit,  DC  voltmeter  for  phase  shift  determination 
and  AC  coupling  to  the  loclcin  amplifier. 

Measurement  of  the  field  strength  across  the  entire  pulse  width 
requires  variation  in  the  relative  delay  between  the  arrival  time  of  the 
pump  and  probe  pulses  at  the  location  of  interaction  within  the  sampling 
head.  The  variation  was  accomplished  by  changing  the  path  length  of  the 
pump  beam  optical  arm  using  a  180  degree  reflector  mounted  on  a  translation 
stage  (Klinger). 

The  signal  periodicity  and  the  reference  waveform  required  by  the 
lock-in  amplifier  for  bandwidth  narrowing  were  provided  by  a  mechanical  beam 
chopper  (Stanford  Research  SR540),  operating  at  2  KHz,  placed  in  the  pump 
beam.  At  each  time  delay  the  lockin  aunplifier  was  allowed  time  to  settle, 
determined  by  the  selected  time  constant,  t,  and  the  output  was  repeatedly 
sampled  and  data  recorded  and  displayed  on  a  computer.  The  signal  to  noise 
ratio  is  improved  by  longer  settling  times  and  larger  number  of  samples, 
however,  this  leads  to  longer  total  experiment  times  during  which  lauser  and 
thermal  instabilities  of  the  electro-optic  medium  effect  the  detected 
signal.  Hence,  there  is  a  compromise  to  be  made  between  elimination  of 
higher  frequency  noise  and  long  term  drift.  Typically  the  time  constant  was 
chosen  to  be  300  ms. 

Characterization  of  the  relationship  between  an  applied  voltage  and  the 
system  response  a  voltage  of  the  form  V^sin(vvt)  was  applied  to  the  LiTaO^ 
modulator,  and  lockin  voltage  recorded.  A  typicad  response  with  slope  of 
4.71  e-3  is  shown  in  figure  11.  Many  applications  of  the  electro-optic 
sampling  system  demand  high  sensitivity.  Hence  calibration  was  extended  to 
to  lower  applied  voltages.  Figure  12  shows  three  examples  of  the  system 
response  in  which  the  lockin  amplifier  output  was  sampled  a  different  number 
of  times  at  each  applied  voltage.  The  response  is  linear  to  within 
experimental  error  and  the  signal  to  noise  ratio  is  seen  to  improve  as  the 
number  of  samples  is  increased  from  1  to  ZOO. 
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V  Applied  (volts) 


Fig.  11  Linear  response  of  the  eiectro- 
optic  sampler  to  applied  bias. 


L/>eont)'  Cot4>roUon 


Fig  12(a) 


Fig.  12(b)  Fig.  12(c) 

Fig.  12  The  signal-to-noise  ratio  is  enhanced  by  increasing  the  number  of 


samples  measui-ed  at  each  applied  bias. 


The  potential  for  signal  to  noise  enhancement  was  determined  by  repeated 
sampling  of  the  lockin  amplifier.  The  applied  voltage  was  held  constant  and 
five  data  sets  of  n  measurements  were  made,  each  measurement  consisted  of 
sampling  the  lockin  amplifier  m  times.  The  average  values  and  standard 
deviation  are  summarized  in  table  1  and  fig  13. 
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1 

2 

3 

4 

5 

19 

19 

19 

9 

9 

1 

10 

20 

100 

200 

5.23 

2.67 

2.92 

1.57 

1.22 

5.75 

5.75 

5.79 

5.79 

5.73 

Data  Set 

n 

m 

<r  (‘leS) 

n-l 

X  (*163) 

Table  1. 


Fig.  13  The  •▼erage  vaiue  and  sample  standard  deviation  for  1,10.20,100,200 

samples  at  constant  applied  bias. 


The  standard  deviation  and  the  slope  were  used  to  determine  the  minimiiTn 
resolvable  voltage  (sensitivity)  for  each  data  set. 

no.  of  samples  1  10  20  100  200 

sensitivity  (mV)  22.4  11.5  12.5  6.5  6.7 


By  recording  200  samples  at  each  applied  voltage  the  error  on  the  measured 
signal  is  reduced  to  6.7  mV  and  the  response  is  linear  within  the  error. 
Recalling  that  the  electrode  separation  is  0.7  mm,  the  minimum  resolvable 
field  is  0.095  V/cm  and  the  minimum  resolvable  E^L  product  is  28.5  e-4 
volts. 


For  short  duration  pulse  field  measurements  it  is  necessary  that  the 
pump  and  probe  beam  path  lengths  result  in  the  same  delay  to  within  the 
accuracy  of  the  pulse  duration.  Careful  measurement  of  the  amn  lengths  is 
essential.  The  synchronous  arrival  of  the  pulses  at  the  sampling  head  was 
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checked  by  temporally  replacing  the  LiTaO^  by  a  cubic  beam  splitter. 
Alignment  was  such  that  both  beams  were  incident  on  the  same  position  of  the 
reflecting  face.  The  transmitted  507.  of  the  pump  beam  and  the  reflected  507. 
of  the  probe  beaun  were  colinear  and  directed  onto  a  fast  photodiode 
(risetime=35  ps),  the  output  of  which  was  displayed  on  an  oscilloscope. 
While  monitoring  the  time  difference  between  the  two  pulses  the  necessary 
minor  adjustments  of  the  arm  lengths  were  made. 

3.1  Electro-optic  Sampling  Measurements 

As  a  test  of  the  electro-optic  sampling  system  a  commercial  photodiode 
(Antel  ARS-2)  was  switched  by  the  pump  beam  and  the  emerging  pulse  launched 
onto  the  balanced  stripline.  The  connection  was  via  50  Q,  27  GHz  SMA 
connectors  and  the  line  was  terminated  with  a  Tektronix  oscilloscope  S4 

sampling  head  (see  fig  14).  Fig.  15  shows  electro-optically  determined  field 
pulses  at  f our  x  positions  in  the  crystal.  The  pulse  is  launched  into  the 

crystal  at  x  =  0  and  measured  at  x  =  1.0,  4.0  6.0  and  8.5  mm.  A  lockin  time 

constant  of  300  ms  was  used  and  at  each  position  200  samples  were  taken.  An 
approximate  measure  of  the  resulting  sensitivity  is  obtained  noting  that  the 
signal  to  noise  ratio  in  the  figure  is  approximately  1  in  60  and  multiplying 
this  value  by  the  150  mV  amplitude  of  the  signed  measured  on  the 

oscilloscope.  The  resolution  is  then  approximately  2.5  mV  which  is  in 
agreement  within  a  factor  of  three  with  the  sensitivity  determiiiad  abor*  in 
section  3.0. 


Fig.  14  A  pulse  launched  from  a  photodiode  is  sampled  by  the  probe  beam  at 
position  X. 
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The  change  in  pulse  shape  as  it  propagates  on  the  line  is  the  result  of 
dispersion  and  reflection.  A  full-wave  analysis  (see  for  example  Carin  and 
Webb,  1989)  using  a  spectral  domain  formulation  was  made  for  a  square  pulse 
of  similar  gross  features  to  that  generated  by  the  photodiode.  The  field 
strength  was  determined  as  a  function  of  time  at  positions  x  =  0  and  x  =  10 
excluding.  Fig  17(a),  aind  including.  Fig  17  (b),  the  effects  of  dispersion. 
The  effects  of  reflection  and  subsequently  dispersion  are  seen  to  be 
significant  in  the  development  of  the  waveform.  These  phenomena  can  be 
expected  to  play  a  role  in  the  evolution  of  the  pulse  propagating  in  the 
balanced  stripline. 

The  characteristics  of  the  electro-optic  sampling  system  are  a  spatial 
resolution  of  less  than  10  |im,  a  sensitivity  of  less  than  7  mV,  and  time 
resolution  of  approximately  7  ps.  The  pulse  widths  demonstrated  are 

approximately  120  ps  FWHM  but  it  is  at  widths  of  <  20  ps  that  the  technique 
offers  measurement  opportunities  not  available  by  other  means.  This  time 
resolution  is  limited  by  the  diiration  of  the  laser  pulse,  not  the 

electro-optic  effect  or  detection  electronics,  hence  the  temporad  resolution 
could  readily  be  improved  by  provision  of  shorter  duration  laser  pulses. 

The  device  has  been  used  to  map  electric  field  pulses  at  severed 
positions  in  a  microstrip  transmission  line  as  a  function  of  time.  The 
measured  pulse  shapes  closely  resemble  those  calculated  using  an  approximate 
model  and  this  application  demonstrates  that  the  technique  could  be  extended 
to  determine  reflection  and  transmission  coefficients  of  pulses  at 

interfaces,  however,  the  purpose  of  this  work  is  to  provide  a  means  of 
studying  transient  field  dynamics  in  high  power  switches.  Hence,  the  next 
step  is  to  incorporate  short  time-scale  semiconductor  switches  with  the 
system. 

The  next  section  outlines  the  phenomena  experienced  by  carriers  at 

metal-semiconductor  contacts,  an  inherent  part  of  any  switch  geometry,  and 
initial  steps  taken  to  fabricate  metal-semiconductor  contacts. 
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g.  IS  Pulse  shapes  as  a  function  of  time  at  (a)  two  x  positions  and  (b) 
four  X  positions. 


(a)  (b) 


Fig.  16  Theoretically  modeled  pulse  shapes  at  x  =  0  and  x  =  10  mm  for  a 
square  pulse  input  both  without  (a)  and  (b)  with  dispersion  included  in  the 
modei. 


4.0  PHOTOCONDUCTIVE  SWITCHES. 

The  switching  times  of  semiconductor  materials  are  determined  by 
carrier  production  and  recombination.  Materials  such  as  GaAs  and 
silicon-on-sapphire  are  known  to  have  recombination  times  of  the  order  of 
several  picoseconds.  Such  materials  appropriately  biased  are  thus  capable  of 
producing  field  pulses  a  few  picoseconds  in  duration.  Hence,  they  provide  a 
more  stringent  test  of  the  electro-optic  sampling  system  or  alternatively 
electro-optic  saunpling  serves  as  a  means  of  investigating  the  switching 
characteristics.  In  designing  a  suitable  switch,  consideration  must  be  given 
to  bulk  and  metal -semiconductor  interface  phenomena  as  the  significance  of 
the  role  of  each  to  current  limitation  varies  with  material,  geometry  and 
applied  bias.  As  an  initial  step  in  the  development  of  switches  for  the 
study  of  carrier  transport  this  section  outlines  the  phenomena  occurring  at 
metal-semiconductor  interfaces  and  describes  some  preliminary  results  of 
f  abrication. 

4.1  Metal  Semiconductor  Contacts 

Metal  semiconductor  contacts  are  inherent  in  photoactivated 
semiconductor  switches.  Schottky,  Stormer  and  Waibel  (1931)  showed  that  the 
potential  drop  in  a  metal-semiconductor-metal  system  occurs  predominantly 
across  the  contact  region.  This  implies  that  there  is  some  form  of  bemrier 
to  conduction  at  the  metal-semiconductor  interface.  Braun  (1874) 
demonstrated  that  the  current  flows  more  easily  in  one  direction  than  the 
other  at  the  contact,  i.e.  the  contact  has  rectifying  properties.  These 
properties  were  exploited  extensively  in  radio  broadcasting  long  before  they 
were  understood.  The  high  resistance  in  one  direction  but  not  in  the  other 
was  presumed  due  exclusively  to  a  layer  of  foreign  materials  at  the 
interface  prior  to  work  of  Schottky  (1938,  1939)  and  Davidov  (1939)  who 
showed  that  asymmetric  behavior  could  arise  as  the  result  of  the 
accumulation  of  space  charge.  It  was  also  realized  that  equilibrium  contact 
conditions  at  every  contact  between  different  materials  involves  charge 
exchange  and  although  quantum  mechanics  predicts  non-linear  behavior  at  any 
lattice  discontinuity  space  charge  effects  dominate. 

4.2  Space  Charge  at  the  Interface. 

Space  charges  are  always  present  at  interfaces  and  this  can  be 
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representatively  demonstrated  by  considering  two  materials  in  contact  and 

having  a  common  current.  Ohm’s  law  states  that  J=<rE  and  if  we  assume  that 

(T  *<r  then  E  ^E  in  the  bulk  sections.  The  non  zero  divergence  at  the 
12  12 

interface  and  Poisson’s  equation  \7.E  =  p/c  leads  to  the  conclusion  ♦’■it  the 
space  charge  density,  p,  at  the  interface  is  non  zero.  Ihis  situation 
results  to  a  double  layer  of  space  charge  at  the  interface  as  illustrated  in 
Fig.  17. 

a,  »o, 

+ _ *  _ - 

— ^  1  hr 
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Fig.  17  Potential  contour  at  any  current-carrying  contact  between  two 
materials  of  different  carrier  concentration  (conductivity).  Note  the 
existence  of  space-charge  double- layer  at  interface. 

The  flux  term  of  the  continuity  continuity  equation,  J  =  e(npE  +  D.Vn), 
describes  the  general  current  density.  For  the  special  case  of  no  gradients 
in  the  carrier  density,  n,  this  reduces  to  Ohm’s  law,  J  =  «rE.  The  diffusion 
current  density,  eD.7n,  is  a  departure  from  Ohm’s  law  and  arises  because  of 
the  space  charge  density.  This  argument  leads  to  the  conclusion  that  third 
phases  are  not  required  for  a  non-linear  resistance  behavior  in  the  contact 
but  can  be  solely  the  result  of  space  charge  at  the  contact. 

4.3  Formation  of  the  Schottky  Barrier. 

Barriers  arise  in  two  distinctly  different  ways.  Firstly  through  the 
difference  in  thermionic  work  f unctions  of  the  two  materials  and  secondly, 
from  the  action  of  surface  states.  Consider  first,  the  difference  in  the 
thermionic  work  fimctions  of  materials  for  the  case  ^  ^  ,  an  important 

1  Rl 

practical  case.  When  the  metal  and  semiconductor  are  separated  the  Fermi 
level  of  the  N  doped  semiconductor  is  greater  than  that  of  the  metal  as  in 
Fig.  2a.  Here  is  the  electron  affinity  of  the  semiconductor.  On  bringing 
the  metal  and  semiconductor  into  thermal  equilibrium  by,  for  example 
connecting  them  with  a  wire  (fig.  18b),  the  Fermi  levels  align  nd  there  is 
a  net  transfer  of  electrons  from  the  semiconductor  to  the  metad.  These 
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electrons  au'e  provided  by  the  donor  atoms  and  as  a  result  of  the  charge 
tramsfer  there  exists  both  a  potential  difference,  V^,  between  the  vacuum 
level  in  the  semiconductor  and  that  of  the  metal  and  a  depletion  layer  of 
width  w  extending  from  the  face  of  the  semiconductor.  The  number  density  of 
donor  atoms  is  orders  of  magnitude  lower  than  the  f  ree  electron 
concentration  in  the  metal,  hence  the  depletion  layer  is  much  thicker  than 
the  excess  chairge  layer  in  the  metal.  As  the  separation  is  decreased  (fig. 
18c)  the  magnitude  of  the  E  field  would  tend  to  increase  without  limit. This, 
however,  is  not  physical  and  is  averted  by  the  further  transfer  of  electrons 
from  the  semiconductor  to  the  metal  until  when  in  contact.  Fig.  18d,  the 
voltage  V  is  zero.  The  potential  S  (  =  d  -  y  )  is  the  barrier  potential, 

l  ns  m  s 

the  energy  an  electron  at  the  Fermi  level  of  the  metal  must  gain  to  be 
raised  over  the  top  of  the  barrier  and  into  the  semiconductor.  The  voltage 
associated  with  the  energy  between  the  bottom  of  the  conduction  band  in 
the  semiconductor  bulk  and  the  top  of  the  barrier  is  known  as  the  diffusion 
potential. 
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Fig.  18  Formation  of  a  barrier  between  a  metal  and  a  semiconductor  (a) 
neutral  and  isolated,  (b)  electrically  connected,  (c)  separated  by  a  narrow 
gap,  (d)  in  perfect  contact,  o  denotes  electron  in  conduction  band;  + 
denotes  donor  ion. 

The  classical  shape  of  the  Schottky  barrier  (1938)  is  parabolic.  This 
can  be  shown  by  assuming  the  semiconductor  is  homogeneous  up  to  the 
boundary,  that  the  uncompensated  donors  give  rise  to  a  uniform  space  charge 
density  in  the  depletion  region  that  is  continuous  rather  than  discrete.  By 
twice  integrating  Poisson’s  equation  for  the  depletion  region  the  energy  of 
the  bottom  of  the  conduction  band  is  determined.  For  N-type  material  of 
permittivity  c,  d\/dx^  =*  -p/e,  where  p  is  the  space  charge  density.  Letting 
E  »  -eV,  be  the  electron  energy,  Poisson’s  equation  can  be  rewritten  as 
d^E/dx^  «  -p/c.  Assuming  donors  in  the  depletion  region  are  completely 
ionized  and  ignoring  both  the  minority  and  majority  carriers,  p  »  eN  ,  in 
the  barrier  and  p  -  0  outside  the  barrier,  we  introduce  a  space  charge 
discontinuity  that  is  wholly  artificial  but  is  a  zero  order  description  of  a 
metal-semiconductor  barrier. 
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d^Z/dx^  =  /e  in  the  barrier, 
d 

On  integration, 

dE  e^Nd  X I  *  =  e^d  (x-X  ) 

dx  e  X  e 
0 

where  X  is  the  barrier  width, 

0 

defined  as  the  value  of  x  for  which 


tne  field  is  zero. 

After  a  second  integration 
2-. 

E(x)  =  - — ^  (x^/2-X  x)  +  const, 
c  0 

With  the  Fermi  level  regarded  as  the 


reference  level, 

E(X  )  =  - — -''o  +  const  =  ^  . 

0  G  n 

•»  const.  »  ^  + 


e^dX^ 

■  ■  0. 
c 


So  E(x)-^  = 

n 


(x-X 

0 


Fig.  19  (a)  Space  charfe  danslty,  (b) 


field,  and  (c)  electron  energy  at  contact. 


e  NdA 

When  X  «  0,  we  have  E(0)-^  «  -= —  o.  This  is  the  barrier  height  seen  from 

n 

eNdX^ 

the  semiconductor  side  and  the  diffusion  potential  is  ^  o. 

The  band  picture,  used  to  explain  bulk  phenomena,  of  a  semiconductor 
arises  because  of  lattice  periodicity.  At  the  surface  of  a  semiconductor  the 
periodicity  is  interrupted  and  the  usual  band  picture  does  not  appiy  at  the 
surface.  As  a  result,  surface  states  exist  and  act  to  trap  electrons  in 
N-type  materials.  The  electrons  are  captured  from  the  donors  in  the  bulk 
before  contact  and  the  donors  which  are  no  longer  compensated  extend  into 
the  semiconductor  a  distance  X^.  The  barrier  thickness  depends  on  the  donor 
concentration  and  the  barrier  height,  thus  the  barrier  thickness  depends  on 
the  energetic  position  and  the  number  of  surface  states. 

The  barrier  discussed  so  far  is  an  abstraction.  It  presupposes  contact 
between  perfectly  clean  materials.  This  in  general  is  not  the  case  of 
experience,  where  each  surface  may  be  covered  with  an  impurity  oxide  layer 
or  the  metal  and  semiconductor  can  interact  with  one  another  to  form  a  third 
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phase.  The  interaction  may  be  inert  or  a  chemical  reaction  may  take  place. 
Often  the  barrier  is  too  thin  (3-20  A)  or  disordered  to  be  studied  by 
standard  band  theory. 


Fig.  20  Metal-semiconductor  contact  with  surface  states. 
The  e  are  four  important  cases  of  third  phases. 


Non-metal 

Metal 

Comment 

(a) 

near  insulator 

physically  adsorbed 

true  work  function 

Schottky  barrier. 

(b) 

covalent 

(eg.  Ge,  Si) 

physically  adsorbed 

work  function  &  surface 

state  dependent. 

(c) 

covalent 

forms  a  compound 

with  non-metaJ 

e.g.  silicide  contacts 

(d) 

covalent 

thin  oxide  film 

e.g.  chemically  etched 

silicon. 

To  a  first  approximation  the  most  important  feature  of  a  contact,  the 
barrier  height  ^  ,  is  a  constant  of  the  two  materials.  For  highly  ionic 

D 

materials,  e.g.  most  group  II- VI  compound  sendconductors,  the  barrier  height 
is  a  strongly  dependent  on  the  work  function  of  the  metal.  These  materials 
have  a  low  density  of  surface  states  and  therefore  the  Fermi  level  is 
unpinned.  In  contrast  group  IV  emd  III-IV  compounds  are  highly  covalent  and 
it  is  an  experimental  fact  that  the  barrier  height  is  for  N-type 


30 


materials  and  for  P-type  materials  where  E  is  the  bandgap.  To  explain 
this  Bardeen  (1947)  proposed  that  the  high  density  of  surface  states  pinned 
the  Fermi  level.  Fig.  21  shows  the  bandgap  and  the  Fermi  level  with  respect 
to  the  bottom  of  the  conduction  band  for  several  semiconductor  materials. 


IPQ$tfifiM  oi  the  twiHhirds  nik  for  barrio  bcighi  pmoiiii 
ai  the  pnetaf-unueondnctor  iaiofact  for  Ao  contacts  on  varioia 
covatem  snacondnctorv  Tlie  tocaiion  of  the  Femi  level  retaiive 
to  the  conduction  band  al  the  interface  is  ploncd  vs 

eocffy  pp  £;  (freo  M.  (4}). 
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Fig.  21  Barrier  energies  for  gold  on  several  n-type  covalent  semiconductors; 
the  straight  line  is  where  the  barrier  energy  is  2/3  of  the  band  gap  (from 
C.A.  Mead.  1966] 

4.4  Current  Transport  Processes  in  a  Forward  Biased  Schottky  Diode. 

Assuming  the  barrier  height  has  been  established,  there  are  four 
mechanisms  via  which  current  can  pass  across  the  interface  independent  of 
barrier  height.  The  current  transport  in  metal-semiconductor  contacts  is 
mainly  due  to  majority  carriers,  in  contrast  to  p-n  junctions  in  which  most 
current  transport  is  by  minority  carriers.  The  mechanisms,  illustrated  in 
Fig  22,  are  (a)  emission  of  electrons  from  the  semiconductor  over  the  top  of 
the  barrier  into  the  metal,  (b)  quantum-mechanical  tuimeling  through  the 
barrier,  (c)  recombination  in  the  space  chamge  region,  (d)  recombination  in 
the  neutral  region  (’hole  injection’). 

4.5  Emission  over  the  Barrier. 

Before  electrons  are  emitted  over  the  barrier  they  must  first  be 
transported  from  the  interior  to  the  interface,  traverse  the  depletion 
region  subject  to  the  mechanisms  of  drift  and  diffusion  and  on  arriving  at 
the  interface  their  emission  into  the  metal  is  controlled  by  the  number  of 
Bloch  states  in  the  metal  communicating  with  states  in  the  semiconductor. 
These  processes  are  effectively  in  series  and  the  current  is  determined  by 
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(a)  Misstoa  of  oitctrons  fnm  tho  Maieon^tor  ovor  tho  top  of  tho 
borrior  into  tho  Mtil» 

(b)  quantta^Mchonical  tunAolllAg  through  tho  barrier. 

(c)  rocoabiaatiM  in  tho  apaco^chargo  rtglMi. 

(4)  roooablnation  In  tho  noucral  region  (*hola  injection'). 


Fig.  22  Transport  processes  in  forward  biased  Schottky  barrier  (from 
Rhoderick  and  Williams,  1988,  p89) 

the  iargest  impediment.  We  assume  that  the  electrons  traverse  to  the 
interf ace  region  the  other  two  processes  have  been  considered  as  current 
limiting  by  Wagner  (1931)  and  Schottky  and  Spenke  (1939)  (drift  and 
diffusion)  and  Bethe  (1942)  (Bloch  states).  The  difference  between  these  two 
models  is  illustrated  by  the  behavior  of  the  quasi-Fermi  level  for  electrons 
in  the  semiconductor.  (The  quasi-Fermi  level  is  a  hypothetical  energy 
introduced  to  describe  the  behavior  of  the  electron  in  non-equilibrium 
conditions.  It  predicts  the  number  of  electrons  in  the  conduction  band  when 
equilibrium  theory  is  assumed  and  the  quaisi-Fermi  level  is  used  in  the 
Fermi-Dirac  distribution  function.) 

In  the  diffusion  theory  the  concentration  of  electrons  in  the 
semiconductor  adjacent  to  the  interface  is  not  altered  by  the  applied  bias 
and  the  quasi-Fermi  level  in  the  semiconductor  is  equal  to  the  Fermi  level 
in  the  metal  (see  Fig.  23).  The  main  obstacle  to  the  flow  of  current  is 
drift  and  diffusion  in  the  depletion  region.  Under  the  assumptions  of 
thermionic  emission  of  Bethe,  electrons  emitted  from  the  semiconductor  are 
not  in  thermal  equilibrium  with  the  metal  conduction  electrons.  The  energy 
exceeds  the  Fermi  level  by  the  barrier  height.  The  bottleneck  here  is  due  to 
the  process  of  emission  of  electrons  into  the  metal. 
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Fig.  23  Electron  quasi-Fermi  level  in  a  forward-biased  Schottky  barrier 

.  according  to  diffusion  theory,  -  according  to  thermionic  emission 

theory.  The  broken  circle  shows  the  energy  distributim  of  electrons  which 
make  their  last  collision  at  a  distance  I  from  the  interface. 


The  current  density  prediction  of  the  diffusion  model  is 

The  current  density  for  thermionic  emission  Is 

J  .  AV«cp[-^j|exp(^)  -  1}  . 

and  that  predicted  by  the  work  of  Crowell  and  Sze(1966)  which  combines  the 


two; 


These  all  have  the  same  dependence  on  4  and  V.  Here,  Vr  is  the 

0 

recombination  velocity,  V  is  the  effective  diffusion  velocity  and  M  ,  is 

4  C 

the  effective  density  of  states  in  the  conduction  band. 


4.6  The  Schottky  Effect 

The  Schottky  effect  is  the  image-force-induced  lowering  of  the 
potential  energy  barrier  for  charge  carrier  emission  when  an  electric  field 
is  applied.  Consider  a  metal-vacuum  system  first.  The  minimum  energy 
required  for  an  eiectron  to  escape  to  the  vacuum  from  the  Fermi  level  in  a 


metal  is  q^  (approx.  2-6  eV).  Where  q^  is  sensitive  to  metal  cleanliness. 

fn  m 


When  an  electron  is  at  a  distance  x  from  the  metal,  a  positive  charge 
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is  induced  on  the  metal  surface.  The  force  of  attraction  between  the 
electron  and  the  induced  positive  charge  is  equivalent  to  force  that  would 
exist  between  an  electron  and  an  equal  positive  charge  located  at  -x.  The 
positive  charge  is  referred  to  as  the  image  charge.  The  attractive  force, 
called  the  image  force,  is 

2 

F  =  -=3_ 

4ne  (2x)^ 

0 

where  is  the  permittivity  of  free  space. 

Work  done  by  an  electron  in  going  from  infinity  to  x  is 


w  =  Jp 


16itc  X 
00  0 


This  corresponds  to  the  potential  energy  of  an  electron  at  x.  When  an 
external  field  is  applied,  the  total  field  is 


PE(x) 


+  qEx  (eV) 


16kc  X 
0 


Fig.  24  is  the  energy  band  diagram  between  the  metal  surface  and  the  vacuum. 


Fig.  24  Elnergy-band  diagram  between  metal  surface  and  the  vacuum.  The  metal 
work  function  Is  qd  .  The  effective  work  function  or  barrier  is  lowered  when 

tn 

the  electric  field  is  applied  to  the  surface.  The  lowering  is  due  to  the 
combined  effects  of  the  field  and  the  image  force.  (After  Sze,  1981) 

The  Schottky  barrier  lowering  (  image  force  lowering)  L4  and  the 
location  of  the  lowering  x  are  given  by  the  condition  dPE(x)/dx  =  0. 
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So 


and 


i.e. 


dPE(x)  ^  -q^  _ 

dx  I6ne 

0  X 


2 


+  qE 


Hence,  A0  =  2x  E. 

m 

These  results  can  also  be  applied  to  a  semiconductor-metal  jimction. 
However,  the  field  should  be  replaced  by  the  maximum  field  at  the  interface 
and  the  free  space  permittivity  replaced  by  the  permittivity  of  the 
semiconductor,  then  =»  (qEy4itc^)  .  The  value  of  may  be  different  from 
the  semiconductor  static  permittivity.  For  GaiAs  the  image  force  permittivity 
for  0<E*'^^<400  (V/cm)^'^^  is  expected  to  be  similar  to  the  static  value. 
Because  of  the  larger  values  of  in  the  metal-semiconductor  system,  the 

barrier  lowerings  are  smaller  than  those  of  a  metal-vacuum  system.  Although 
the  barrier  lowering  is  small  it  has  a  significant  effect  on  the  current 
transport  processes  in  metal-semiconductor  systems. 


4.7  Tunneling  (or  Field)  Emission  Through  the  Barrier. 

Under  some  circumstances  electrons  with  energies  below  the  top  of  the 
barrier  penetrate  by  quantum  mechanical  tunneling.  The  current  density  is  of 
the  form 


E  is  the  diffusion  potential  of  the  Schottky  barrier  and  is 
00 


where  m*  is  the  effective  mass  of  an  electron. 

The  current  density  is  proportional  to  (transition  probability)  x 
(occupied  probability  in  semiconductor)  x  (unoccupied  probability)  in  metal. 
If  the  temperature  is  raised  the  electrons  are  excited  to  higher  energy 
levels  and  see  a  thinner  and  lower  barrier,  consequently  the  current  density 
increases,  this  higher  temperature  phenomena  is  known  as  thermionic  field 
emission  and  is  illustrated  in  Fig.  25. 


35 


Fig.  25  Field  and  Thermionic  field  emission  under  forward  bias.  The  diagram 
ref  ers  to  a  degenerately  doped  semiconductor  f  or  which  ^  is  negative. 
(Padovini  and  Stratton.  1966) 


4.8  'Ohmic*  Contacts  and  Field  Emission 

Field  emission  is  important  to  production  of  'Ohmic*  contacts.  Ohmic 
contacts  do  not  necessarily  obey  Ohm's  Law  but  are  of  low  resistance 
compared  to  that  of  the  bulk.  They  are  usually  Schottky  barriers  on  very 
highly  doped  materials.  The  expression  for  the  contact  resistance,  Rc,  is  in 
general  a  complicated  expression,  Yu  (1970),  but  the  doping  dependence  is 


for  field  emission 


where  E 

00 


for  thermionic  field  emission 

for  thermionic  emission. 


Hence  for  the  field  and  thermionic  and  field  emission  the  current  increases 


at  higher  doping  levels. 


4.9  Fabrication  Techniques. 

Ohmic  contacts  can  be  fabricated  by  doping  the  semiconductor  in  the 
contact  region.  Fabrication  requires  the  deposition  of  a  metal,  or  mixture 
of  metals,  and  a  dopant  followed  by  an  annealing  process.  The  deposition  is 
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by  one  of  a  variety  of  techniques  including  evaporation,  sputtering, 
pressure  contacts,  soldering  or  plating.  We  chose  evaporation.  In  producing 
a  contact  one  is  faced  by  a  several  obstacles.  The  wetting  action  of 
metallic  compounds  on  GaAs  is  low.  The  contact  materials  ’ball  up’  on 
cmnealing  forming  small  islands  on  the  substrate  and  resulting  in  a  poor 
contact.  When  using  a  mixture  of  metals  it  is  necessary  to  evaporate  to 
completion  to  maintain  the  correct  stoichiometry.  The  steps  in  the  procedure 
are 

•  etch  the  GaAs  substrate  to  remove  oxide  layers  on  the  surface 

•  evaporate  and  deposit  Au-Ge, 

•  evaporate  and  deposit  Ni. 

•  anneal  to  promote  interdiffusion, 

•  evaporate  and  deposit  a  metal  conducting  strip. 

Tunneling  electron  microscopy  indicates  that  the  presence  of  a  NiAs(Ge)/GaAs 
interface  (Fig.  26)  is  important  to  lower  the  contact  resistance  and  that 
carrier  transportation  via  tunneling  in  Ge  doped  GaAs  could  be  most 
important  Ogawa  (1980)  has  concluded  that  high  reliability  is  realized  when 
the  Ni  layer  is  thinner  than  the  AuGc  and  that  a  rapid  anneal  to  500  C 
improves  the  depth  of  alloying  and  the  contact  character  (Fig.  27). 


MMMIMC  TtMPtMATunC  T  (T  I 


Fig  26  Schematic  illustration  of 
the  sequence  of  alloying  reactions. 
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Fig.  27  Alloying  depth  vs.  annealing 
for  two  different  times  at  maximum 
temperature. 


Back  to  back  Schottky  diodes  were  fabricated  amd  a  chopped,  focused  laser 
beam  was  directed  at  one  contact  and  then  the  other  (Fig.  28)  the  resulting 
output  current  was  oppositely  directed  when  changing  from  one  contact  to  the 
other  as  would  be  expected  (Fig.  29).  This  indicates  that  Schottky  barriers 
were  indeed  fabricated. 


Fig.  28  Geometry  of  gold  electrodes 
on  GaAs  substrate  illuminated  by 
chopped  laser  beam. 


Fig.  29  Current  flowing  under 
illumination  of  one  contact 
(top)  and  the  second  (bottom). 


Under  the  assumptions  that  light  changes  the  barrier  by 

1.  changing  the  boundary  carrier  concentration  on  the  bulk  side, 

2.  photoexcitation  of  traps  within  the  barrier  region  and  subsequent  removal 
of  carriers  by  the  high  field, 

3.  the  adjustment  of  the  trap  occupancies  within  the  space  charge  region  to 
the  new  majority  and  minority  carrier  concentrations, 

Bardeen  (1950)  and  Branbury  (1953)  obtained  an  expression  for  the  open 


circuit  voltage 


V 

oc 


+  1 


where  I  is  the  saturation  current,  I  is  the  current  due  to  incident 

•  1 
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Photovoltaic  Voltage 


photons  and  n  is  the  non-ideality  factor  of  the  diode.  Figure  30  shows  the 
voltage  measured  by  a  1M£J  input  impedance  oscilloscope  across  the  contact  as 
a  function  of  the  incident  beam  power.  The  slope  of  the  line  for  sample  2  is 
0.030  and  kT/q  =  0.026  at  300  K.  The  slope  of  the  curves  varies,  i.e.  n 
varies,  n  may  take  values  between  1  and  2.  n  =1  when  diffusion  currents 
dominate  and  2  when  recombination  currents  dominate.  The  curve  for  sample  6, 
the  sample  which  was  prepared  by  evaporation  to  completion,  rapidly  annealed 
and  on  LEG  GaAs  has  the  lowest  photovoltaic  voltage  and  lowest  value  of  n 
indicating  that  it  is  the  most  ’Ohmic'  of  the  contacts  and  that  diffusion 
currents  are  more  important. 


Fig.  30  Photovoltaic  response  of  gold-gallium  arsenide  contacts  to  laser 
illumination. 

4.10  Chromium  Doped  GaAs  Switch  Response. 

Fiirther  experimentation  with  methods  of  depositing  Au/Ge/Ni  contacts  to 
Cr:GaAs  provided  improved  contacts  as  defined  by  the  voltage  switched  when 
photoactivated  and  the  time  response  of  the  switching  circxiit.  These 
contacts  made  feasible  the  measurement  of  the  switch  response  to  short 
duration  illiunination  by  electro-optic  sampling.  A  detailed  description  of 
the  technique  for  contact  deposition  is  provided  in  appendix  B. 

To  improve  the  transmission  line  characteristics  over  that  described 
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earlier  where  the  fast  photodiode  was  external  to  the  LiTaO^  electro-optic 
sampling  head,  Cr:GaAs  was  incorporated  in  the  into  the  transmission  line  as 
shown  in  Fig.  31.  The  electrodes  continued  from  the  LiTaO^  across  the 
Cr:GaAs  except  for  a  small  (25  pm)  gap  in  the  top  electrode  1  mm  from  the 
LiTaO  /CrtGaAs  interface.  An  external  bias  of  6  volts  established  a  field 

3 

of  2.4  KV/cm  across  the  gap  which  was  illuminated  from  above  by  180  fs 
FWHM,  76  MHz,  150  mW,  620  nm  laser  pulse  train. 


Fig.  31  Cr:GaAs  switch  incorporated  with  LiTaO^  balanced  stripline. 

This  geometry  reduced  both  the  number  of  discontinuities  and  the 
distance  through  which  the  switching  pulse  propagated  before  being  sampled. 
The  LiTaO^  was  reduced  to  0.35  mm  resulting  in  a  reduced  half  wave  voltage 
of  2.95  KV.  The  measured  data  points  are  compared  in  figure  32  with  a 
theoretically  determined  waveform.  This  waveform  assumes  that  the  voltage 
pulse  decreases  exponentially  with  a  (fitted)  lifetime  of  20  ps  as  expected 
due  to  charge  carrier  recombination  and  that  it  is  reflected  from  the 
Cr:GaAs/LiTaO  with  a  coefficient  of  .2886,  determined  from  the  mismatch  in 

3 

the  transmission  line  impedance.  The  time  between  successive  reflections  is 
59  ps  as  determined  from  the  position  of  the  sampling  beam  from  the  edge  of 
the  gap  and  the  velocities  of  the  pulse  in  the  two  media.  Although  the  peaks 
align  well  in  time  the  amplitudes  and  the  shape  vary  considerably  after  the 
first  pass  of  the  pulse.  Possible  reasons  for  these  departures,  yet  to  be 
confirmed,  are  that  the  calculated  reflection  coefficient  is  not  accurate. 
This  may  arise  from  the  use  of  Hammerstad’s  expressions  for  the 
determination  of  the  characteristicimpedance  in  materials  of  relative 
dielectric  constant  greater  than  16.  A  value  for  which  their  accuracy  is  not 
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known.  Another  obvious  shortcoming  of  the  calculated  curves  are  that  they 
madce  no  inclusion  of  dispersion  which  is  significant  for  pulses  containing 
high  frequency  components. 


Fig.  32  Electro-opticallj  nKMunored  Cr:CaAs  twitch  retponae  to  ItO  fs  FWHM 

laser  pulse. 


5.0  Conclusion 

An  electro-optic  sampling  system  has  been  demonstrated  for  measurement 
of  electric  fields  having  variations  on  ultra-short  timescales.  The 
technique  has  been  used  to  measure  the  response  of  both  commercially 
available  photodiodes  and  the  photoresponse  of  biased  Cr:GaAs.  The  field 
switched  in  the  Cr:GaAs  is  resolvable  because  of  the  lowered  impedance  of 
the  metal-semiconductor  contacts.  The  contacts  ar  the  result  of 
establishing  a  deposition  and  annealing  procedure  for  AuGe.'N’i  on  Cr:GaAs.  In 
the  absence  of  such  a  procedure  the  photoresponse  would  not  be  measurable. 

The  technique  will  now  be  applied  to  the  detailed  investigation  of 

elc''tron  dynamics  in  high  power  semiconductor  switches,  including  the 

dependence  of  switching  behavior  on  electric  and  magnetic  fields,  and  the 

position  of  illumination.  These  investigations  will  be  directed  towards  the 

understanding  of  the  persistant  photoconductive  phenomena  known  as  Tockon’. 
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APPENDIX  A 

Pockels  Cell  Analysis  via  Jones  Calculus. 

The  state  of  polarization  of  a  beam  transmitted  through  numerous 
optical  elements  can  be  determined  using  Jones  calculus.  Assume  the  light  to 
be  normally  incident  on  the  systems  elements,  then  it  is  possible  to 
represent  the  effect  of  any  element  on  the  polarization  or  intensity  as  a 
linear  operator  acting  on  the  electric  vector  or  the  light  wave.  The 
operator  is  represented  in  the  form  of  a  2  x  2  matrix.  In  general  the  effect 
of  an  element  with  different  refractive  indices  and  absorption  coefficients 
along  the  principle  axes  of  the  element  can  be  treated.  However,  we  assume 
absorption  is  zero. 

Assume  the  y  direction  is  out  of  the  page  and  we  are  looking  down  the 

y  direction  from  a  location  at  larger  y  than  any  element  and  represent  the 

field  in  two  components. 

_  _  1(0  _  _  _  ^1(0 
E  =  E  e  x  and  E  =  E^  e 

X  Ox  z  Oz 


Which  we  abbreviate  to 

E  =  E 


and  E 


E. 


X  Ox’  ~  X  Oz 

The  state  of  polarization  of  the  output  of  an  optical  element  is  related 
to  the  input  via  a  matrix  known  as  a  Jorws  matrix.  For  example,  consider  the 
change  in  the  polarization  state  brought  about  by  a  lineair  polarizer 
orientated  at  angle  0  with  respect  to  the  x-y  coordinate  system  of  Fig.Al  . 


Fig.  A1  E  components  separated  by  the  polarizer  rotated  to  angle  0. 


2. 

COS  0 

sin0  COS0' 

fE  ] 

X 

COS0  sin0 

•  2  . 
sin  0 

The  Pockel’s  cell  Jones  matrix 
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representation  is 


E’' 

1 

'cos3 

sinp' 

1  0  ' 

'l-2a  1 

fE  ' 

x 

o  ‘T 
10  e 

X 

E’ 

“  2 

-sin^ 

cos^ 

1  l+2a 

V  / 

■^zJ 

(1)  (11)  (111) 


where  (i)  represents  the  rotation  function  of  the  second  polarizer  at  angle 
P  with  respect  to  the  z  axis,  (ii)  the  total  phase  shift,  and  (iii)  the 
polarization  matrix  of  the  first  polarizer  oriented  at  angle  a  +  n/4  radians 
from  the  z  axis,  where  a  is  small.  After  multiplication 


E’  ^  =  il(E^+E^)-4a(E^cos^^-E^sin^^)+sin2^.cos^(E^+E^)-2asin2^cos^(E^-E^)] 

x4xzx  z  xz  xz 

E’  ^  »  i[(E*+E^)-4a(E^sin^^-E^cos^^)-sin2^.cos^(E^+E^)+2asin2^cos^(E^-E^)]. 

Z4X2X  Z  xz  xz 

'2  *2 

These  expressions  indicate  that  E^  and  E^  are  complicated  functions  of  a 
and  ^  and  that  it  is  of  importance  to  set  a  »  0  and  ^  »  x/4.  W.  n  a  -  0 

E’^=  i{(EWKl+sin2^.cosr)l 

X  4  x  z 

E'^=  i((E^+E^)(l-sin2^.cos^)] 

z  4  X  z 

and  when  p  =  -it/4 

E*^=  iUE^+E^Kl-cosDl  =  ilE^+E^)sinV/2 

X  4  X  z  2  X  z 

E’^»  il(E*+E*)(l+cosr)I  -  i(E^+E^)cos^/2 

z  4  X  z  2  X  z 

hence  E'^-E’^*  i(E*-*-E^)lsinV/2  -  cosV/2l. 

X  z  2  X  z 

E’^-E’^»  i(E^+E^)(l-2cos*(r/2)J 

X  z  2  X  z 

Taking  the  difference  between  the  two  signals  leads  to  a  doubling  in  the 
measured  voltage  swing  for  a  given  change  in  phase  shift.  For  small  f  both 
E^  and  E'  have  a  linear  dependence  on  the  introduced  phase  shift,  Fig.  Al. 

2  2 

E’  amd  E’  are  directed  to  separate  diodes,  and  subsequently  to  separate 

X  z 

inputs  of  the  lockin  amplifier  differential  input.  Using  the  diffentiad 
input  in  this  manner  both  increases  the  signal  measured  and  enables  an 
increased  sensitivity  measurement  to  be  made. 
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Appendix  B 


Fabrication  of  GaAs  Switches 


Materials 

1.  A  doped  substrate  with  one  side  polished 


5mm  (0.196  in.)  L  x  0.135  mm  (0.053  in.)  w  lapped  to  035  mm  (0.014 
in.)  thick 


impurity: 

concentration 

resistivity: 

mobility: 

thickness: 

Surface  direction: 
Estatic 


Cr 

SxlO^Vcm^ 

>10^  ncm 
4000  cm^/ vs 
originally  0.020  inches 
2®  oflF(lOO) 

=  12.8 


Effp  —  11.0 


2.  Li  Ta  O3  Orientation  z  parallel  to  0.7  mm  dimension. 

10  mm  L  X  13  mm  (0.059)  W  x  0.7  mm  (0.027  in)  thick 
10  mm  6x13  mm  (0.054  in.)  w 

X  0.7  mm  (0.027  in.)  +  0.07  mm  W  Faces  with  optical  polish,  flat, 
parallel 

13  mm  W  Faces  hand  lapped  and  polished 


3.  Mica  Backing 

OJ0038*  thick 


.08"  2inm — H 
.04" 


~  2.5mm  .1" 


_5minJ 
.197" 


-1- 


4.  All  12%  Ge 

0.0015  in.  thick  pieces  cut  to  weigh  out  exact  amount  for  500-700  A 
thick  film  when  evaporated  to  completion  from  Tungsten  (W)  boat. 

5.  Nickel  wire  99.97%  purity  0.010  in.  diameter 
10  in.  wound  on  Tungsten  filament 

6.  Tra-Con  Epoxy 

Tra-bond  BA-2115 

MAJOR  STEPS 

I.  Deposition  of  Au  12%  Ge  (500-600  A)  -I-  Ni  (125-150  A) 

II.  High  temperature  (510-540®C)  Alloy 

III.  Deposit  aluminum  (8(XX)A)  electrodes  on  backs  of  GaAs  -t-  LiTaOa 

IV.  Epoxy  GaAs  end  to  LiTaOa  on  mica  backing 

V.  Deposit  aluminum  (8000A)  on  epoxy  junction  and  over  mica  and  ends 

VI  Deposit  aluminum  (SOOOA)  electrodes  on  gap  side  of  assembly  with  gap  masked 

Dff. 

The  following  is  a  detailed  description  of  the  above  procedures. 

Deposit  Au  12  %  Ge  to  completion  —  6(X)A  thick  -  on  GaAs 
I  A.  Vacuum  System  -  Set-up 

1.  Difiiision-pumped  vacuum  system 

2.  Tungsten  (W)  "dimple”  boat  source  with  0.0315  gms  Au  12%  Ge  pieces 

3 

3.  4—  in.  throw  -  substrate  placed  directly  over  source 

4 

3 

4.  Glow-discharge  ring  l--in.  below  substrate 

4 

5.  Quartz  crystal  oscillater  thickness  monitor  next  to  substrate 

6.  Microscope  slide  plate  for  thickness  measurement  on  profilometer  (Dektak)  IIA 
next  to  substrate, 
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Microscope  slide  next  to  substrate.  Used  for  thickness  measurement  on 
(DEKTAK IIA)  after  completion  of  run. 

7.  Resistance  heater  over  substrate 


IB.  Clean  and  Set-up 
1.  Clean  substrate 

Ultrasonic  clean  2  mins,  each  in  trichloroeth^ne;  acetone;  meth^ 
alcohol,  D  J.  H2O  ,  blow  dry  vith  nitrogen.  Chemical  surface  etdi: 

2  mins.  (3)  H3PO4:(l)H2O2:(50)  H2O 
5  mins.  20ccNH4OH:7ccH2O2:97ccH2O, 

Deionized  H2O  rinse  and  dry  nitrogen  blow  dry 


2.  Immediately  load  in  mask-substrate  holder 

The  holder  centers  the  GaAs  0.054”  polished  face  over  a  0.015  in.  wide  slot.  A 
0.001  in.  diameter  nichrome  wire  stretched  across  the  slot  creates  a  gap  in  the 
AuGe  Film  contact  strip  about  0.040  in.  from  the  end  of  the  GaAs 


Substrate  holder  mounted  directly  over  Tungsten  boat  and  Au  Ge  source  to 
prevent  any  shadowing  of  wire  masking.  Resistance  heater  sits  over  holder. 
Thermocouple  on  holder  to  monitor  temperature. 

IC  Thin  Film  Deposition  Procedure 

Pump  down  to  <  2xlO'^T 
Heat  up  to  -  150®C 
Cool  down  to  -  80“C 

Glow-Discharge:  3  mins.  40mA  0.85kv  43mT  N2  96®C 
Pump  down  to  <2xlO‘^T 


-3- 


Deposition: 

Deposit  to  completion-to  assure  proper  alloy  of  Au-Ge-until  thickness 
monitor  shows  no  change 

~  4  kHz  change  During  Deposition  -  -  5(X)-600Ath‘'  Typically  80-90 
wig  130  sup  {o}  c  .br 

Dry  N2  up  to  atmospheric  pressure. 


DEPOSIT  TOP  NICKEL  FILM  -  USA 

Replace  W  boat  +  Au  Ge  with  W  filament  -I-  Ni 

GaAs  over  Ni  filament 

Without  removing  from  holder,  place 

Proceed  as  in  lA  step  C  up  to  deposition 

Deposition:  400  Hz  change  on  monitor  (shutter  controls) 

slow  DEPOSITION  ~  30  sec. 

Cool  to  <70®C(DryN2  ATM)- 
EXAMINATION  OF  FILM  - 


1.  Check  width  of  Gap  ai^  "  Quality”  on  digital  measuring  microscope.  If  no 
shadowing  and  clean  edges:  ~  25n  wide  then 


2.  Take  picture  with  filars  at  edges  of  Gap, 

3.  Check  thickness  of  both  AuGe  and  Ni  on  DEKTAK  -  print  out  AuGe: 
~  500-600A  Ni'.-  125-150A 


II.  Alkqr-In 


Contact  deposition  for  OHMIC  contact  to  GaAs  "Rapid  Thermal  Process":  Using 
Tantalum  Boat  for  heating  and  ooolinf 


Indent  for 
GaAs 

0.005”  Tantahim 
Ta  cover  with  hcries  over  top  to  confine  heat  but  let  in  nitioaM 


re- 


-A- 


Nitrogen  fed  in  through  needle  valve  in  base  plate  Weight  on  bell  jar  to  allow  positive 
pressure  of  nitrogen 

A.  Pump  dovm  •  N2  flush 

Roughing  pump  and  N2  decrease  pressure  to  25  mT. 

N2  increase  pressure  to  ~  1  mm  of  Hg  Several  times  to  flush  out  other  gases 
Raise  to  atmospheric  pressure  and  then  POSITIVE  PRESSURE  with  KN2  Gas 
flow  through  a  needle  valve 

B.  Heat  up  in  N2  ATM 


Increase  the  voltage  across  the  tantalum  boat  in  steps. 


5V  20  secs. 

lOV  40in.  temperature  increased  to  340°C 

13V  55in.  564®C 
-  12V  5  mins.  30  secs.  540-509®C 
(stayed  mostly  -  517®C ) 


-OFF- 

OV  1  min.  63°C 
Summary:  55  secsat  temperature  increased  to  364°C 

5  min.  509-540®c  ( -2ymins.  at  517®c)|  positive  press  N2  atmosphere 

1  min.  temperature  decreased  to  63”C 

C.  Inspect: 


Measure  Gap  -  Take  Picture 

Gap  should  be  same  width  as  before  alloying  - 

Should  be  no  large  globs!! 

Very  small  Xjeometric"  shapes  appear  to  be  preferentially  oriented 
par^l  to  Gap,  across  width  of  GaAs 
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III.  Deposit  8OOOA  Aluminum  Electrodes  on  Back  of  GaAs  and  LiTa03 

A.  Vacuum  System 

Diffusion  pumped  system, 

Sputter  gun  with  aluminum  target. 

Throw  -  Horizontal  Gun  to  Substrates 
Glow  -  Discharge  -  Ring 

-  No  Heater  - 


B.  Set'Up 


-  No  -  Cleaning 

Back  side  of  GaAs  and  0.054  in.  wide  hand-polished  face  of  LiTa03  over  0.015  in.  W 
slot  in  holder 


H — Li  T  aO  GaAs  1*0— 


Holder  clamped  on  opposite  gun 

C.  Deposition 

Pump  to  <  3.0x10^7 


Gap  Up 


Mica  Underneath 


Mca 


Glow  Discharge;  3  mins.  25ma.  85kV  -  90mt  Argon  Pressure 
Deposition  (Pre-Sputter  5  min.) 

FILM  DEPOSITION:  8  min.  630V-770  mA  3.0  mT  Argon  Pressure 
8  mm.  630V  -  770  ma  3.0mT  Argon 
cool  -30' 

Argon  increase  pressure  to  1000  mT  and  then  introduce  air  and  increase  pressure  to 
ATM 


IV.  EPOXY  GaAs  TO  UTaO}  ON  MICA  BACKING 

A  special  jig  is  used  to  line  up  the  GaAs  and  LiTa03  end  to  end  on  a  mica  backing 
piece.  The  mica  piece  is  notched  at  the  junction  so  a  thick  aluminum  film  may  be 
deposited  after  gluing  to  bridge  the  gap  in  the  back  side  electrodes  deposited  in  Step 
III.  Teflon  ')[)ushers”  from  the  side,  top  and  end  assure  flatness  and  alignment. 
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1 .  Very  light  Him  of  silicon  release  agent  wiped  on  jig 

2.  Small  dab  of  Tra-Con  Tra-Bond-21  IS  put  on  either  side  of  notch  in  mica 

3.  Mica  back  placed  in  jig 

4.  Very  small  dab  of  epoxy  put  on  end  of  LiTa03 

5.  LiTaOs  now  placed  in  jig,  epoxy  end  centered  over  notch  in  mica 

6.  GaAs,  with  gap-end  nearest  junction,  now  carefully  placed  in  jig 

7.  Side,  top  and  end  '^ushers*  tightened  in  place 

8.  Inspection  --  Jig  is  open  underneath  at  junction  so  it  may  be  viewed  under 
microscope 

9.  Cure  epoxy:  over  night  or  4  hrs.  at  65®C 

V.  DEPOSIT  ALUMINUMSOOOATHICK  ON  EPOXY  JUNCTION  -SAME 
PROCEDURE  AS  STEP  III 

1 .  Clean  mica  and  ends  carefully  with  trichloroethylene  on  swab. 

2.  Place  epoxied  assembly  over  O.OIS  in.  wide  slot  so  aluminum  is  deposited  on 
back  side  -  at  junction,  on  mica  and  eiqwsed  ends. 

3.  Deposition  as  in  III  C. 

VI.  DEPOSIT  ALUMINUM  ELECTRODES  SOOOA  THICK  ON  SIDE  GAP. 

1.  The  epoxied  assembly  is  now  placed  on  the  O.OIS  in.  wide  slotted  jig  so  that  the 
gap  is  centered  over  a  0.02S  wide  bridge  which  masks  off  the  gap.  The 
electrodes  are  deposited  on  either  side. 

2.  Deposition  as  in  III  C. 
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Abstract  -  A  non-equilibrium  multi-vaiiey  transport  model  for 
carriers  in  compound  semiconductor  devices  has  been  developed. 
This  macroscopic  transport  model  provides  an  efficienc  scheme  for 
device  modeling,  and  can  overcome  the  difBcuity  in  evaluating  the 
relaxation  times  in  multi-valley  conservation  equations  without  a 
priori  assumption  of  the  displaced-Maxwellian  distribution.  This 
model  has  been  successfully  applied  to  modeling  of  electron 
transport  in  GaAs  subjected  to  a  rapidly  time-varying  6eld. 


where//*, 7  is  the  distribution  in  the  ith  valley,  N is  the  number  of 
valleys,  *,■  is  the  wave  vector,  £  is  the  electric  field,  and  h  is  the 
Plank  constant  divided  by  2w.  The  first  term  of  the  RHS  repre¬ 
sents  the  intravalley  scattering,  and  the  intervalley  scattering  pro¬ 
cesses  out  of  and  into  the  ith  valley  are  described  by  the  second 
and  third  terms,  respectively.  The  carrier  density,  energy  and 
momentum  conservation  equations  in  each  valley  can  be  obtained 
from  Eq.  (1)  by  averaging  over  the  distnbution,  and  are  written  as 

ch  /V  nr 

♦  V-(ii,u,)  -  53  -  53  n^M^,  (2) 


Most  of  transport  models  that  have  been  used  (l]-[6]  to  ac¬ 
count  for  non-equilibrium  phenomena  in  compound  semiconductor 
devices  is  based  on  the  single-valley  conservation  equations.  It  has 
been  assumed  that  the  energy  and  momentum  relaxation  times  (r, 
and  r„,  respectively)  only  depend  on  mean  energy  (theoretically, 
on  all  moments  (7]  of  the  distribution)  and  is  described  by  the 
phenomenological  equations  [3|,[4|.  In  compound  semiconductor 
devices,  carrier  behavior  is  described  by  the  distribution  in  each 
valley.  The  distributions  in  different  valleys  obey  the  Boltzmann 
transport  equations  (BTEs)  that  are  coupled  by  intervalley 
scattenng.  Consequently,  the  macroscopic  description  is  given  by 
a  set  of  multi-valley  conservation  equations.  These  equations  in 
different  valleys  are  thus  coupled  by  intervalley  traiuitions. 
However,  the  relaxation  times  in  multi-valley  conservation  equa¬ 
tions  can  not  be  determined  from  the  phenomenological  equations 
[3], [4].  The  difficulty  in  determining  these  times  has  led  to  the  fact 
that  the  use  of  conservation  equations  has  assumed  only  a 
single  valley  where  the  carrier  density,  mean  energy  and  average 
momentum  are  averaged  over  all  valleys. 

In  order  to  determine  the  relaxation  times  and  to  account 
for  intervalley  transfer  effects,  an  assumption  has  been  made  in 
most  papets  [8]-(12]  using  multi-valley  conservation  equations  that 
the  distribution  function  in  each  valley  is  displaced  Maxwellian. 
These  distributions  are  then  used  to  evaluate  the  relaxation  times. 
However,  this  assumption  is  rather  severe  since  the  distribution  in 
compound  semiconductors  is  appreciably  different  from  a  displaced 
Maxwellian. 

In  this  paper,  a  transport  model  based  on  the  concept  of 
the  Macro-ldnedc  distiibiitioa  (MKD)  function  (7],(13],(14]  (whose 
time  scales  are  characteriaed  by  momencs  of  the  distribution)  has 
been  developed.  The  description  provides  an  efficient  scheme  for 
device  modeling,  and  can  account  for  intervalley  transfer  effects 
without  a  priori  assumption  of  a  displaced  Maxwellian  distribution. 


The  BTE  in  each  valley  is  given  by 


-  *  ^  V. 


^1  ,53  f^i  . 


V(«,P,a,) 


qEn,  -  V(a,*,r^)  ~  -ifl  -  53  -ill. 


V(a,«,)  +  V  [n,a,(a,  ♦  kp'J]  -  qS  n^m, 

.  -A  - «.) .  ii  .  £  v., 

y-uw  y-uw 

where  kg  is  the  Boltzmann  constant,  and  is  the  electron 
temperature  related  to  carrier  energy  ^  by 

m’  is  the  macroscopic  effective  mass  of  carnets,  and  n,,  and 
Uj  are  carrier  density,  average  momentum,  and  average  velocity  in 
the  ith  valley,  respectively.  The  averaged  intervalley  carrier  transi¬ 
tion  rates  (v^)  are  defined  as 

(5) 

where  is  the  transition  rate  from  the  ith  to  jth  valley.  The  intra- 
valley  and  intervalley  momentum  exchange  times  and  r„y, 
respectively)  are  given  by 

(6a) 


'tiL. .  - 


and  the  intravalley  and  intervalley  energy  exchange  times  (r^, 
r^  and  r^)  are  defined  as 


-  «.) 


)/(»')  -  Oa) 
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-r 

'•Vi 


(e.  *  Ae^,  T  hu^^)f^(kj) 


and  the  information  in  the  time  scale  smaller  than  r„  Is  filtered. 
/m«  i^^n  be  obtained  from  Eq.  (9a)  in  terms  of  or  r„  and  the 
relaxation  times  or  rates  in  each  valley  can  be  evaluated  from  the 
integrals  given  in  Eqs.(5)-(7)  using/^,,,.  These  rates  or  times  are, 
therefore,  tabulated  for  solving  the  multi-valley  conservation  equa¬ 
tions  given  in  Eqs.  (2)-(4). 

(7c) 

m.  Results  and  Conclusion 


with 

The  energy  separation  between  the  jth  and  ith  valleys  is  A€j  =E,j  - 
E,„  where  £  (E^)  is  the  energy  g^p  between  the  jth  (ith)  valley 
and  valence  band,  and  huji  is  the  intervalley  phonon  energy. 

In  the  previous  work,  [71,[13],[14)  the  concept  of  the  MKDs 
has  been  introduced  to  describe  carrier  behavior  in  various  time 
scales.  Depending  on  the  device  scale,  a  MKD  in  each  valley  cor¬ 
responding  to  the  scale  of  interest  is  chosen  to  evaluate  the  relaxa¬ 
tion  times.  The  coarseness  of  the  macroscopic  description,  there¬ 
fore,  depends  on  the  choice  of  the  MKD  fi^.  The  choice  is  based 
on  the  fact  that,  generally,  in  semiconductors,  where  r 

is  the  carrier  density  relaxation  time.  For  example,  if  the  scale  of 
interest  is  in  the  order  of  r,,  the  MKD  is  chosen  to  be/j^,(k)  = 
/*,,(kn,e)  whose  time  scale  is  described  by  r  and  r,;  thus 
contains  only  information  to  the  scale  greater  than  or  as  small  as 
T,.  Similarly,  in  the  scale  smaller  than  r,  (therefore  nea^  t„),  a 
higher-resolution  MKD  can  be  taken  as  f^Jk) ~fM„(k,n, « ,  f).  In 
this  paper,  only  the  encrgy-tirae-scale  MUD  f^,  is  discussed.  The 
momentum-tirae-scale  MKD  f^m  studied  in  the  future. 

In  the  energy  time  scale,  the  time  scale  of  the  carrier  distri¬ 
bution  in  the  ith  valley  can  be  changed  [7],(14)  in  the  form: 

/,(*,r,r)  -  (r.r),e,(r,r)]  (8) 


The  developed  multi-valley  transport  model  has  been 
applied  to  the  study  of  electron  transport  in  GaAs  subjected  to  a 
step  electric  field  of  finite  rise  time.  The  velocity  response  to  the 
step  change  in  Beld  is  shown  in  Fig.  1  where  included  are  results 
obtained  from  the  developed  multi-valley  transport  model,  Monte 
Carlo  simulations  using  a  multi-valley  band  model,  and  the  single¬ 
valley  approximation.  The  average  velocity,  shown  in  Fig.  1, 
obtained  from  the  single-valley  conservation  equations  is  very 
similar  to  the  previous  results  given  in  [5];  Le.  the  single-valley 
approximation  leads  to  faster  relaxation  than  the  Monte  Carlo 
simulations.  This  is  due  to  the  fact  that  at  the  applied  field 
electron  behavior  strongly  depends  on  the  intervalley  scattering. 
The  single-valley  approach  can  not  properly  account  for  slower 
averaged  intervalley  transitions,  and  thus  results  in  faster  response. 
However,  as  can  be  seen  in  Fig.  1,  the  multi-valley  transport  model 
can  take  into  account  the  intervalley  transfer,  and  therefore 
provides  a  better  description  than  the  single-valley  approach. 


Fig.  1.  Average  velocity  of  electrons  in  n-type  GaAs 
subjected  to  a  step  change  in  field. 


That  is,  in  the  r,  scale,  the  space-time  dependence  of/,  is  taken 
to  result  from  changes  in  n„  and  t,  Using  Eq.  (8)  in  the  BTE, 
given  in  Eq.  (1),  to  change  the  scale  of  evolution  of  the  distribu¬ 
tion,  it  leads  to  a  "macro-kinetic"  equation  for  the  MDK  in  the 
scale  of  r,.  The  equation  for  the  r, -scale  MKD  has  been  found 
to  be,  for  field  along  z  direction 


with 

U(  is  the  c-vi.-ier  velocity.  Eq.  (9a)  represents  a  distribution  that 
obe-^  a  steady-state  BTTE  with  a  time-varying  "equivalent"  field  of 
the  ith  valle  /,  £j^  as  the  driving  field,  and  E,^  is  given  by  Eq.  (9b) 
in  terms  of  t,.  This  MKD  thus  responds  instantaneously  to  Ej^ 
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